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Using a sample of high-redshift lensed quasars from the CASTLES project 
with observed-frame ultraviolet or optical and near-infrared spectra, we have 
searched for possible biases between supermassive black hole (BH) mass esti- 
mates based on the C IV, Ha and H/3 broad emission lines. Our sample is based 
upon that of Greene, Peng & Ludwig, expanded with new near-IR spectroscopic 
observations, consistently analyzed high S/N optical spectra, and consistent con- 
tinuum luminosity estimates at 5100A. We find that BH mass estimates based 
on the FWHM of Civ show a systematic offset with respect to those obtained 
from the line dispersion, 07, of the same emission line, but not with those ob- 
tained from the FWHM of Ha and H/3. The magnitude of the offset depends 
on the treatment of the He II and Fe II emission blended with C IV, but there is 
little scatter for any fixed measurement prescription. While we otherwise find 
no systematic offsets between Civ and Balmer line mass estimates, we do find 
that the residuals between them are strongly correlated with the ratio of the 
UV and optical continuum luminosities. This means that much of the disper- 
sion in previous comparisons of C IV and H/3 BH mass estimates are due to the 
continuum luminosities rather than any properties of the lines. Removing this 
dependency reduces the scatter between the UV- and optical-based BH mass 
estimates by a factor of approximately 2, from roughly 0.35 to 0.18 dex. The 
dispersion is smallest when comparing the Civ 07 mass estimate, after remov- 
ing the offset from the FWHM estimates, and either Balmer line mass estimate. 
The correlation with the continuum slope is likely due to a combination of red- 
dening, host contamination and object-dependent SED shapes. When we add 
additional heterogeneous measurements from the literature, the results are un- 
changed. Moreover, in a trial observation of a remaining outlier, the origin of the 
deviation is clearly due to unrecognized absorption in a low S/N spectrum. This 
not only highlights the importance of the quality of the observations, but also 
raises the question if whether cases like this one are common in the literature, 
further biasing comparisons between C IV and other broad emission lines. 

Subject headings: gravitational lensing — galaxies: active — quasars: emission 
lines 



1. Introduction 

It is thought that every massive galaxy has a supermassive black hole (BH) at its 
center, and some physical properties of the BH appear to be tightly correlated with those 
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of the galaxy. In particular, the mass of the central BH correlates we l l wit 



Tremaine et al. 2002: Gultekinet al. 2009 



Ferrarese fc Merritt 2000; Gcbhardt et al. 



i the luminosity 



of the spheroidal component of the host (see, e.g.. iMarconi fe Hunt 12003 
and with its velocity dispersion (see, e.g. 



Graham! 12007) 



2000 



Graham et al.l 201ll ) . Both of these properties of 



galaxies have physical scales a few orders of magnitude larger than the sphere of influence of 
the BH, so mechanisms linking their properties are not immediately apparent. Theoretical 
models try to account for the correlation through co-evolution of the galaxy and its BH, in 
which accretion induced by galaxy mergers regulates the BH's growth, and feedback from 
the accretion r egulates the growth of the galaxy by quenching star formation and removing 



cold gas (e.g., Granato et al.l 12004; iHopkins et al.ll2005l 120061 120081 : ISomerville et al.l 12008 



Shankar et al.ll2009l ). However, the existence of these correlations does not necessarily imply 
co-evolutionary mechanisms, as some aut hors argue that they can be a simple cons equence 
of mergers and the central limit theorem (jPengl 120071 . l2010l : iJahnke &: Maccioll2010l ). 



Direct measurements of BH masses in inactive galaxies are only possi ble for a small num- 
ber of nearby objects becau se it is necessary, or at least desirable (see, e.g.. lMerritt fc Ferrarese 
200ll : iGiiltekin et al.ll2009l ). to resolve the BH's sphere of influence in order to determine the 
BH mass from the kinematics of the stars and gas closest to it. Galaxies with active nuclei 
(AGNs) offer a completely different means of estimating BH masses at any distance. In 
particular, Type 1 AGNs show bright broad emission lines in their spectra produced by gas 
in the broad line region (BLR), which is close to the central black hole but outside the hot 
accretion disk. The large line-widths are thought to arise from the Doppler broadening due 
to the orbital velocity of the gas around the BH. Thus, measuring the mass of the central 
BH from the width of the broad lines is possible if the distance of the BLR from the BH is 
known. 

T his distance can be directly measured wit h the reverberation mapping (RM) tech- 
nique falandford & McKee! Il982l : IPetersonl Il993l ). This technique works by measuring the 
light travel time between the continuum and the broad-line emitting regions, which is de- 
rived from the time lag between changes in their respective luminosities. Unfortunately the 
times cale over which appreciable variability is observed in AGNs increases with BH mas s 
(e.g., IVanden Berk et al.l [20o3 : Iwilhite et aJbood : iKelly et al.lbood : iMacLeod etaJboioh . 



making it difficult (i.e., more time inte nsive) to apply RM to t he luminous QSOs that possess 
the most massive BHs. For example, IMacLeod et al.l ( 120101 ) find that for a typical quasar 
with Mbh = 10 s M (typical magnitude of Mj ~ —23 mag), the rest-frame timescale, At, 
required to reach an r.m.s. variability amplitude of 0.1 mag is approximately 45 days, while 
for a quasar with Mbh — 10 9 M Q (Mj rj —25.5 mag), At is approximately 125 days. This is 
further complicated by the time dilation due to the higher redshift of these rare objects. It 
has been shown, however, that the distance from the BH to the BLR correlates well with the 
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continuum lum inosity of the AGN (see, e.g-. lKaspi et aI.ll2000L 120051 : Bentz et alJl200d 120091 ; 
Zu et al.ll2010l ). Given this correlation, BH masses can be estimated for distant broad-line 
quasars for which the RM technique is not reasonably applicable. Masses estimated in this 
way are usually referred to as single epoch (SE) BH mass estimates. 

Because it is generally easier to obtain optical rather than UV or IR spectra, SE BH 
masses are typically estimated from the H/3 and Ha broad emission lines and the con- 
tinuum luminosity at 5 100 A at low redshifts (z < 0.7). The overlap with RM targets 



has allowed for very accurate calibration of these SE m ass estimators (jCollin et al. 2006 



Vestergaard & Peterson 2006. IVP06l ; iMcGill et al. 20081 ). At high redshift, however, these 
emission lines are shifted into the IR, and most mass estimates are then based upon the UV 
MgllA2798 and CivA1549 broad emission lines and the continuum luminosities at 3000A 
for Mgll and 1450A or 1350A for Civ. Unlike the Balmer emission lines, these UV lines 
lack large l ocal calibration samp l es be cause of the difficulty of obtaining UV-based RM mea- 



imp . 

surements. IQnken fc Kollmeierl (j2008l ) have argued that the Mg II line can provide accurate 
mass estimates, but that there is a small, but significant, dependence on the Eddington 
ratio of the AGN. C IV, on the other hand, is not thought to have this bias, and IVP06 
have calibrated a Civ-based mass estimator based on local RM AGNs using space-based 
UV spectra. However, there are still concerns about whether the Civ velocity widths are 
attributed solely to gravity or if there are bulk flows due to winds of ejected material, and 
the impact of these effects on the accuracy of C iv-based BH mass estimates is still debated. 
For example, Civ is typically slightly displaced in wav elength (usual l y blueshifted) with 
respect to the rest o f the quasar emission lines (see, e.g., Gaskell 1982 ; Tvtler fc Fan! 19921; 
Richards et al. 2002 ). and frequently shows broad absorption features (e.g., I Wevmann et al. 



1981 ) and strong line asymmetries correlat ed with quasar properties (e.g., 
Richards et aDl2002[ iLeighlv k Moore! l2004h - 



Wilkes 



1984 



The simplest approach to test the reliability of C IV mass estimat es is to systematicall y 
compare them to Balmer line estimates for the same sources (see, e.g.. iDietrich et al.ll2009l ). 
High redshift lensed quasars are some of the best targets for such tests. Generally, the prob- 
lem is that the high redshift makes it easy to observe the C IV line, but the better calibrated 
Ha and H/3 lines lie in the near-IR, where it is difficult to observe them. Magnification 
increases the apparent brightness of the lensed quasars, and also, because their observed 
brightness is not uniquely determined by their intrinsic luminosity and distance, it helps 
to mitigate any Eddington biases in the sample or, in other words, it makes objects in the 
sample unlikely to be preferentially brighter than average for their BH mass. 



In a recent work, Greene, Peng & Ludwig (2010. iGPLlOhl presented near-IR spectral 
observations for a sub-sample of lensed quasars from the CfA- Arizona Space Telescope LEns 
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Survey (CASTLES) of gravitational lenses ( Falco et al. 
had been estimated in a previous work by iPeng et al. 



20011) w hose Cr y or Mgll BH masses 



( 120061 ). iGPLlOl measured, whenever 



possible, the full width at half maximum (FWHM) of the H/3 and Ha emission lines of these 
objects and found no systematic biases between BH masses estimated from these lines and 
those estimated from C IV. Their sample, however, did not cover a large enough range in 
BH mass to decide whether there was a mass dependent slope to the rela tion between the 
masses. This comparison also suffered from the fact that IPeng et al.l (120061 ) lacked access to 
the original UV/optical spectra for many targets and frequently had to rely on the printed 
spectra in published papers to measure line widths. 



In this work we start from the sample of lGPLlOl and attempt to improve on both of these 
issues. First, we add Balmer line based BH mass estimates for the lens SDSS1 138+0314 and 
make revised estimates based on new, higher S/N, spectra of HS0810+2554 and SBS0909+532. 
We obtained near-IR observations for SDSS1138+0314 and HS0810+2554 using the newly 
commissioned Large Binoc ular Telescope (LBT) NIR Spectrograph Utility Camera and Inte- 



gral Field Unit (LUCIFER; ISeifert et al.ll2003l: lAgeorges et a 



we use the UV through IR observations of Mediavilla et al. ( 201ol ). Second, we made consis 



20101), while for SBS0909+532 



te nt C IV BH mass es timates from high S/N spectra using the original observations analyzed 
by IPeng et al.l (120061 ). other published or unpublished spectra, or new spectra for all targets 
in the sample. Finally, we obtained continuum luminosities at 5100A for all objects in 
the sample in a consistent manner. This allows us to i nclude t he lenses SDSS0246-0825, 
HS0810+2554 and Q2237+030, which were excluded by IGPLlOl . With these additions we 
expand the sample of IGPLlOl with both C IV and Balmer lines mass estimates from 7 to 
12 quasars and the mass range covered by approximately 0.5 dex. In £J2] we describe the 
sample of gravitationally lensed quasars we use in this study as well as our observations. 
In £J3] we describe the methods we use to measure emission line velocity widths and their 
uncertainties, the continuum luminosities of the quasars and the SE BH masses. In §4] we 
compare the different mass estimates we have derived and determine the possible biases we 
measure between them while in £J5] we expand our results using a heterogeneous sample of 
measurements from other studies. In £0 we summarize the conclusions. In an appendix we 
discuss individual objects in detail. We use a standard ACDM cosmology with Q M = 0.3, 
fl\ = 0.7 and H = 73 km s" 1 Mpc -1 throughout the paper. 



2. The Sample of Lensed QSOs 



We selected 12 lensed quasars from the CASTLES survey with high quality UV/optical, 
typically ground-based, spectra of C IV and either published near-IR spectra of the Balmer 
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lines or IR magnitudes bright enough to obtain such spectra. The targets are listed in Table 
[TJ All 12 objects have been observed by CASTLES with HSTm the U(F555W), /(F814W) 
and H (F160W) bands, except for B1422+231, which was not observed in I. 



We start from the sample of iGPLlOl . who observed most of these lensed quasars in the 
near-IR with the Triplespec spectrograph at the Apache Point Observatory. The wavelength 
range of these spectra is 0.95-2.46 fim with R = 3500, and either the B. 8 or Ha (or both) 
emission line is observable in one of the atmospheric windows. Although IGPLlOl considered 
objects with a large span of redshift and reddening, we limit our sample to objects with 
sufficiently high redshift and small enough reddening for C IV emission to be observable in 
ground based UV/optical spectra0. IGPLIO presented FWHM velocity width measurements 
for all the objects in their sample but did not present BH mass estimates for three of them. 
For these three lensed QSOs, SDSS0246-0825, HS0810+2554 and Q2237+030, we have mea- 
sured the continuum luminosity and estimated BH masses so we can include them in our 
sample. 

We obtained near-IR spectra in the if and K band for SDSS1138+0314 (Fig. [T]) and in 
the J-band for HS0810+2554 (Fig. [2} with the LBT LUCIFER spectrograph. The first was 
obtained as part of the LUCIFER science demonstration time an d is discussed h ere, whil e 
the second was a target of a separate project to be presented by iMogren et al.l ( lin prep.l ). 
We also analyzed the near-IR J- and if-band observations of SBS0909+532 presented by 
Mediavilla et al.l (J2010J) , shown in Figure [3J 



2.1. LUCIFER Observations of SDSS1138+0314 

We obtained a near-infrared spectrum of SDSS1138+0314 using the new LUCIFER 
instrument at the LBT during its science demonstration time. LUCIFER is a near-infrared 
spectrograph and imager with an overall wavelength range of 0.85 - 2.5 /im. We observed 
SDSS1138+0314 in the longslit mode with the OrderSep filter, a 0"5 slit, the 200.H+K grating 
and the Nl . 8 camera for a total integration time of 840s over 7 dithered exposures during 
the night of UTC 2010-01-04. This configuration gives an effective wavelength range of 1.49 
- 2.4 /im, which includes both the if and if bands, with a resolving power of 1880 at if and 
2570 at K. The slit was oriented to include images A and C of the lensed quasar, as well 



] Peng et al.l |2006h mistakenly quote a C iv BH mass estimate for the lens J1004+1229 that, in fact, 
corresponds to the lens SD SS1004+4112 (C.Y. Peng, private communication). This error was propagated 
into the analysis of IGPLlOl . The quasar in J1004+1229 is highly reddened and it is not possible to see its 
Civ emission in a UV/optical spectrum. 
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as part of the lens galaxy. No emission from the lens galaxy is detected in our data. The 
B9V star HIP 33350 was observed with the same configuration, except for a change in slit 
width from 0'.'5 to 1", and was used to correct the spectrum of SDSS1138+0314 for telluric 
absorption features. The difference in resolution caused by the different slit widths degrades 
our telluric corrections, but has little consequence for measuring the width of broad emission 
lines. We estimated the seeing was ~ 0'/8 during the observations. 

We reduced the d ata using standard IRAF packages in combination with the IDL task 
xtellcor_general of IVacca et al.l ( 120031 ) for the telluric absorption corrections. We per- 
formed a 2-D wavelength calibration on each of the 7 exposures using the sky emission lines 
and built a sky frame by median combining them. The sky frame was then used to remove 
the sky from each exposure before extracting the spectrum. We also did an alternate sky 
subtraction of the spectra using a version of the COSMOS software modified to work on 
LUCIFER data. Thi s software, designed for reduction of sp ectral observations with IMA CS 



(IDressler et al.l 120061 ) and LDSS-3 (upgraded from LDSS-2 , Allington- Smith et al.lll994T ) on 
the Magellan telescopes, follows the procedures of iKelsonl ( 120031 ). It produces an accurate 
model of the sky emission by creating a sub-pixel resolution map of the sky line profiles 
using the full extent of the lines in the spectrum coupled with a model of the optical dis- 
tortions. Both extractions of the spectra yield equivalent results, and both are shown in 
Figure [TJ While in principle we could use the telluric standard to perform an absolute flux 
calibration, it is hard to model the slit losses, especially considering the difference in the slit 
widths. Instead, the flux calibration was performed by convolving the spectrum corrected 
for telluric absorption with the NICMOS F160W filter curve and matching it to the esti- 
mated de-magnified absolute magnitude of the quasar from the CASTLES HST imaging of 
this lens (see §3.31 for details on the lens magnification). The blue edge of the LUCIFER 
SDSS1138+0314 spectrum is somewhat redder than the blue edge of t he F160W b a nd, so 
we extended the observed spectrum using the AGN SED template of lAssef et al.l (120101 ) 
assuming no reddening. Note that the HST NICMOS observations were obtained on UTC 
2003-11-06, approximately 6 years before the LUCIFER observations, so we attempt to cor- 
rect for the intrinsic variability of the quasar. However, this is typically not an important 
correction (see §3.3p . We use the R— band light curves o btained with the SM ARTS 1.3 m 



telescope for a gravitational lens monitoring project (see iMorgan et al.ll2010l ). These data 
show that the quasar intrinsically brightened by 56 ± 17% between UTC 2004-02-03 and 
UTC 2010-01-09. We assume that no significant variability occurred between the HST NIC- 
MOS and the first SMARTS observations and between the LUCIFER and the last SMARTS 
observations. 



From an optical spectrum of SDSS1138+0314, lEigenbrod et al.l (120061 ) estimated a red- 
shift of z — 2.438 for the quasar, while SDSS provides z = 2.4427 ±0.0014. Using the narrow 
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component of Ha and the [OIII] AA4959, 5007 emission lines, we obtained z = 2.4417, con- 
sistent with SDSS. We did not use the [Nil] lines or the narrow component of H/3 as they 
could not be centroided accurately because of blending with the broad Ha and H/3 profiles, 
respectively 



2.2. UV/Optical Spectra 



For most of the IGPLIOI sample, as well as for SDSS1 138+0314, we found suitable high 
S/N optical spectroscopic observations in the literature that the owners kindly made avail- 
able for this study (see Table [T] for the references, where applicable, and Appendix |A] for 
details on each object). When needed, we performed an absolute flux calibration using pho- 
tometry from several different sources, as this was not always required for the science goals 
of the original project. All the UV/optical spectra compiled from the literature are shown 
in Figure HI 

We could not locate suitable optical spectra for HS0810+2554 and FQB1633+3134. 
Both objects were observed by the SDSS spectroscopic survey, but these spectra did not have 
high enough S/N to provide accurate line- width measurements with good continuum sub- 
traction. We obtained new optical spectra of these objects using the MDM observatory 2.4m 
Hiltner telescope with the Boiler & Chivens CCD Spectrograph (CCDS). HS0810+2554 was 
observed on UTC 2010-02-24 with a grating center of 5300A and was flux calibrated using the 
standard star Feige 34. FBQ1633+3134 was observed on UTC 2010-03-21 and UTC 2010- 
03-22 with a grating center of 4700A and was flux calibrated using the standard star Feige 
98. Absolute fluxes were obtained for both objects by performing a cross-calibration between 
SDSS g-band photometry of other o bjects in the field an d g-band photometric observations 



with the RETROCAM instrument (Morgan et al.l 120051 ) obtained on UTC 2010-03-06 and 



UTC 2010-03-22 for HS0810+2554 and FBQ 1633+3134, respectively. The reduced spectra 
are shown in Figure HJ 



3. Models and Measurements 



In this section, we briefly discuss the methods we use to measure the line widths and 
estimate the black hole masses from the optical and near-IR spectra. 



2 http : //www . astronomy . ohio- state . edu/MDM/CCDS/ 
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3.1. Line- width Measurements 



There is no standard prescription for measuring the line-width characterizations of the 
broad C IV emission line in QSOs. While for other emission lines this may not be a significant 
source of uncertainties, there is a shelf-like emission feature redward of C IV that blends with 
the line profile and is created by a combination of broad HeIIA1640, OIII] A166 3, and a 



feature of unknown origin at 1600 A usual ly referred to as the A1600 feature ( Laor et al. 



1994 ; iMarziani et al.lll996l; Fine et al.ll2010l) . While the A 1600 feature is commonly thought 



to correspond to Fe II, Fine et al.l (l2010h argue that this cannot account for all the observed 
flux, yet it is also unlikely that C IV can reach large enough velocities to produce the feature. 
Different prescri ptions for modeling the blended emis sion can have significant effects on line 



width estimates fD ennev et al 



2009J Fine et al.ll2010f ). so it is important to explore how these 



affect our results. iFine et al.l ( )2010l ) explored three different and widely used approaches and 
their effects on the C IV width measurements. The three prescriptions are: (1) to assume that 
the A1600 feature corresponds to C IV emission and therefore remove only the He II A1640 and 
OIII] A 1663 contributions; (2) to assume that the A1600 feature belongs to a different species 
from Civ and so removing its contribution along with that of the other two c omponents 



on th e shelf; and (3) to fit the A 1600 featur e as part of the continuum (see IFine et al. 



2010l . for details on each prescription). While IFine et al.l (120101 ) selects prescription (2) as 
their preferred method, in large part because it produces symme tric C IV profiles, it is hard 
to apply this approach to low S/N data (see IFine et al.l |2010| . for details). Moreover, it 
is not guaranteed to produce more accurate BH masses than the other two prescriptions. 
The simple prescription of (3) produces line-width characterizations that are systematically 
smaller than prescription (2) but with very low dispersion between individual measurements, 
while (1) produces estimates with a larger scatter relative to (2) but without a systematic 
offset. The differences between the prescriptions is smallest for FWHM and largest for the 
line dispersion, <j[. 

Based, in part, on these issues, we considered two different prescriptions for removing 
the continuum and blended emission from the C IV emission line profile. Both prescriptions 
are amenable to large scale automated use. The first p rescription, which we will refer to 
as prescription A, is very similar to that used by |VP06| . where the shelf feature redward of 
Civ is considered part of the Civ line profile, but only the region within ±10,000 km s _1 
of the peak is considered. The continuum is fit by linearly interpolating between the two 
continuum windows in the wavelength ranges 1425-1470 and 1680-1705A. When these con- 
tinuum windows were affected by absorption, we slightly s hifted t hem as detailed in Table |2j 
Our continuum fitting is in principle different from that of IVP06l . who considered 5 different 
continuum windows and then fit a power-law to them, but the differences of the measured 
line-widths are not significant and our approach requires a much smaller wavelength range 
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for the spectra. In this prescription, He II and O ill] emission is not explicitly removed, but 
this has negligible effects due to th e limit on the velocity range, making it analogous to 
prescriptio n (1) of iFine et al.l ( 120101 ). The second prescription, B, is analogous to prescrip- 
tion (3) of IFine et al.l (120101 ). as we fit the A1600 feature as part of the continuum. It only 
differs in that the red continuum region is chosen to match the minimum between C IV and 
the A1600 feature. In general, prescription A will lead to broader estimates of the Civ line 
width than prescription B. 

The observed wavelength continuum windows for each object and prescription are listed 
in Table |2j The C IV emission line flux was then measured above the fit continuum and 
between the emission line wavelength regions listed in Table [2j In addition, for objects that 
showed mild absorption features, bad pixels, and/or significant night sky line residuals, we 
used a low-order polynomial (i.e., first, second or third order depending on the size and 
location of the feature) to interpolate across the feature before measuring the line widths. 
Details for the individual targets are given in Appendix |A] We did not attempt to remove 
any narrow-line emis sion from C IV A 1 549, since this line is typically ver y weak and cannot 
be reliably isolated (j Wills et al.l Il993l although see ISulentic et al.l 120071 ) , and the separate 



VP06 



and references therein for 



lines of the C IV doublet are unresolved in AGN spectra (see 
further discussion). We characterized the line width by both its FWHM and line dispersion 
(ai, the second moment of the line profile). The widths were measured directly from the 
actual or interpolated spect rum (except where no ted below and in Appendix [A]) following 
the procedures described by iPeterson et al.l (120041 ) . 



We also fit the original or interpolated line profiles with a sixth-order Gauss-Hermite 
(GH) polynomial, because making functional fits to emission-line profiles i s a common way 
of mitigating the e ffects of low S/N on line-width measurements (see, e.g.. Woo et al.l 120071 ; 



McGill et al. 



20081. for similar approaches) . The Gauss-Hermite polynomial s we fit utilize the 



norma lization of Ivan der Marel fc Franxl (119931 ) and the functional forms of ICappellari et al. 
(120021 ). We then use a Levenberg-Marquardt least-squares fitting procedure to determine 



the best-fitting coefficients. We measured the widths of these line profile m odels using the 



same software as was used to measure widths directly from the data (see IPeterson et al. 



20041 ). Ultimately we only used the results from the line profile models for PG1115+080 (see 



Appendix |AJ. Instead, these fits were primarily used to determine uncertainties in our width 
measurements as described in £13.21 The continuum and the Gauss-Hermite fits to the C IV 



line profiles are shown in Figure [5] for both prescriptions. In the cases of SDSS1138+0314 
and SBS0909+532, reasonable fits could not be achieved because of the extremely high S/N 
and peculiar shape of these line profiles (a very narrow peak with broad base; see Appendix 
|A|) . No fits are shown for these objects. 
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Both the FWHM and line dispersion, 07, measurements of the CivA1549 emission line 
are listed in Table [3] for all objects in our sample for both prescripti o ns. W e have corrected 
the widths for spectral resolution effects following iPeterson et al.l (120041 ). when possible, 
using the resolutions given in Table |5J Except for PG1115+080, we utilize the line widths 
measured directly from the data (interpolated across gaps where noted) for the subsequent 
black hole mass calculations. For objects with multiple spectra of the individually lensed 
im ages we averaged their line widths. Our Civ A 1549 widths are smaller than those given 
bv lGPLlol for the objects in which we both used the SDSS spectra (Q0142-100, SDSS0246- 
0825, PG1115+080, and H1413+117). The likely origin of the discrepancy is that |GPLlo| fit a 
narrow line component as part o f the C I V profile, which would naturally yield larger FWHM 
values. We note, however, that GPL10 do not use their SDSS line-width measurements to 



estimate BH masses in their analysis, but always use those determined by lPeng et al.l ( 120061 ). 
The lens HE 1 104-1805 i s the only obje ct in the sampl e for w hich we use the same optical 
spectrum as IPeng et al.l ( 120061 ). that of IWisotzki et al.l (119951 ). and we find a FWHM that 
is smaller by 260 km s _1 , compared to our measured uncertainty of 50 km s -1 . Although 
Peng et al.l (120061 ) do not quote errors in their line width measurements, the disagreement 
(~ 5 A in the observed-frame) is likely within their uncertainties. 

Line widths of the H/3 and Ha broad-emission lines are given in Table [31 while the 
continuum and broad line spectral wavelength regions used are given in Table HJ We mea- 
sured them from the near-IR spectra following a similar procedure to the Civ line- widths 
except that (1) the best Gauss-Hermite polynomial fit was used for all line-width measure- 
ments, with the exception of Ha for SBS0909+532, because the S/N of the near-IR data 
was typically too poor to justify measurement directly from the data, (2) blended emission- 
line components were removed from each spectrum before the line width was measured, as 
described in Appendix [AJ and (3) a power-law, instead of a linear, continuum was fit to 
the H/3 spectrum of HS08 10+2554 because it was fit simultaneously with additional blended 
emission-line components over a larger wavelength range. 

For the objects where we lack the nea r-IR spectroscopic observations, we rely on the 
published Ha and H/3 line widths of lGPLlOl These measure ments we re done using somewhat 
different methods than ours. While we consider most of the iGPLlOl FWHM estimates to be 
reliable, there are some that we believe are suspect because (1) they were measured from 
very low S/N spectra, (2) the lines were not fully contained in the wavelength range of the 
spectrum, and/or (3) we do not agree with the narrow-line component models subtracted 
before the line width was measured. In the relevant Figures and Tables, we differentiate 
between the Balmer-line velocity widths we think are reliable (group I, solid symbols) and 
those we believe are affected by any of these issues (group II, open symbols). Individual 
objects can be in both groups because these issues may affect only one of the Balmer lines. We 
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also include in group I the Ha and H/3 line- width measurements from our new IR spectra. The 
decision to split our sample is a conservative choice, and our conclusions are not significantly 
modified when the group II line widths are included. 



Line-width measurements can be affected by sources of error that are difficult to model, 
as they depend not only on the overall S/N ratio, but also on the line profile and the presence 
of sky emission and absorption lines, with the latter being of particular importance in the 
near-IR. We use a Monte Carlo approach to determine the uncertainties in our line-width 
measurements. Using the flux uncertainty per pixel in each spectrum and the best fit Gauss- 
Hermite line profile (with the exception of the optical SDSS1138+0314 and SBS0909+532 
spectra, see Appendix |A]) , we produced 1000 resampled spectra by adding random Gaussian 
deviates based on the error spectrum to the flux in each pixel of the GH model spectrum and 
then re-measured the line width using the methods described in the previous section. For 
the UV/optical spectra from the literature without an error spectrum, we estimated one by 
propagating the measured S/N of a small continuum window near the CivA1549 emission 
line to the overall spectrum. In this case, 5 Fx, the flux error in a pixel of wavelength A with 
flux Fx, is given by 



where A c and F\ c are the average wavelength and flux per unit wavelength of the continuum 
window chosen, and (S/N) c is the signal-to- noise ratio per pixel in the chosen continuum 
window. This equation is constructed by assuming that the only source of error is Pois- 
son fluctuations, and that the number of detected photons is proportional to Fx(hc/ A) -1 , 
where the proportionality constant is empirically determined in the continuum window from 
(S/N) c , A c and Fx c . This approach neglects the sky background and the presence of strong 
absorption or emission sky lines, which is reasonable for the UV/optical spectra. It also 
neglects changes in the instrument sensitivity as a function of wavelength and assumes a 
constant pixel wavelength- width, both of which are reasonable because the continuum S/N 
is measured in close proximity to the emission line of interest. While the parametric fits 
are not exact representations of each line, this still provides a reasonable estimate of the 
fractional uncertainties. 



3.2. 



Line- Width Measurement Uncertainties 
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3.3. Luminosity Measurements 



We estim a ted th e continuum luminosities at 5 100 A by fitting the AGN SED template 
of lAssef et al.l (120 10f ) to the unmagnified quasar magnitudes obtained from the CASTLES 
project HST NICMOS imaging. To correct the observed quasar fluxes for the lens magnifi- 
cation, we modeled each system using the astrometry and lens galaxy photome try from the 



CAST LES iZST WFPC2 and NICMOS observations following the procedures of iLehar et al. 
(120001 ). The image is decomposed into a set of point sources for the quasars, de Vau 



couleurs models for the lens galaxy and, if necessary, a lensed host component, convolved 
with model or empirical P SFs. The resu lting component positions and image fluxes were 
modeled using lensmodel ( iKeetonl 120011 ). The lens was modeled as a singular isothermal 
ellipsoid in an external shear with the ellipsoid's orientation and ellipticity constrained by 
those of the light of the lens galaxy and a weak prior on the external shear. The models 
were not tightly constrained to match the observed fluxes due to systematic errors in im- 
age flux ratios such as source variability and microlensing. Aside from substructure, the 
dominant uncertainty in the magnification s is the radial mass distribution of the lens (see 
Kochanek. Schneider fc Wambsgansd 120041 ). and this is less than a factor of two even if we 
allow the full range of models between a flat rotation curve and a constant Mj L model. Since 
we have extensive evidence tha t lenses have mass distributions corresponding to flat rota- 



tion c urves on these scales (e.g.. lRusin et al.ll2003t : I Jiang fc Kochanekll2007t Koopmans et al. 



20091 ). the model uncertainties are considerably less than this factor, and the uncertainties 
are dominated by the systematic uncertainties in the image fluxes. Table [1] lists the magnifi- 
cations used for each object in the sample. The only object for which a different model w as 
used is Q0957+561, where we used the magnifications determined by iFadely et al.l ( 120101 ) . 



We did not apply reddening corrections other than removing Galactic foreground extinc- 
tion (see below), as the requirement that Civ is observable in the UV/optical severely limits 
the presence of dust absorption, especially at rest-frame 5100A. For all four-image lenses, 
we estimated the true source flux for all images, rejected the highest and lowest estimates 
and averaged the remaining two to limit the effects of microlensing. For two-image lenses 
we simply averaged the two estimates. Table [1] shows the estimated unmagnified H-hand 
magnitude of each quasar. Note that in general we did not apply a correction for var iability. 
Although there is a 5 to 10 year time difference between the CASTLES and the IGPLIO 
Triplespec observations, the typical uncertainty introduced falls well below the systematic 
uncertainties in the SE BH mass estimates. An estimate of the typical variability of a quasar 
can be obtained from mea surements of their structure function. Using the power-law fit of 
Vanden Berk et al.l (120041 ) to the i-band structure function of SDSS quasars, we find that 
the typical quasar would experience a change in magnitude of approximately 0.2 mag for 
a rest-frame time-lag of 1500 days (approximately 10 years in the observer's frame for our 
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lowest redshift quasar). A change of 0.2 magnitudes results in a change to the BH mass 
estimate of 0.04 dex, well below their typical error bar of 0.3 dex, and we would expect the 
H— band variability to be still s maller, as the average variability amplitude de creases with 
increasing wavelength (see, e.g.. lVanden Berk et al.ll2004l ; iMacLeod et al.ll2010l ). 



For SDSS1138+0314, HS0810+2554 and SBS0909+532, we performed an absolute flux 
calibration of the near-IR spectra and measured the 5 100 A continuum luminosity directly. 
The calibration for the first object is discussed in detail in §2.11 For HS08 10+2554 we fit a 
power-law to the continuum of our MDM CCDS spectrum (see §2.21) and extrapolated it to 



rest-frame 5100A. For SBS0909+532 we calibrated the spectrum using the HST NICMOS 
i/-band photometry, as the object did not show significant flux variations between the two 
relevant epochs (J. Munoz, private communication). 

To obtain the rest-frame continuum UV luminosities at 1350A and 1450 A, we flux cal- 
ibrated the spectra whenever it was necessary and measured the flux by fitting a straight 
line to the region between rest-frame 1349A and 1355A for the estimate at 1350A and to 
the region between 1440A and 1460A for the estimate at 1450A. We corrected these lumi- 
nosities for foreground Galactic ex tinctions obtained thr ough the NASA/IPAC Extragalactic 
Databas^fl from the dust maps of Schlegel et al. ( 1998 ). Errors in the continuum luminos- 
ity will be dominated by the uncertainties in the magnification models, which are hard to 
quantify. We assume a conservative error of 20% in each continuum luminosity estimate. 



GPLlOl obtained continuum luminosities at 5100A for their sample of objects by follow- 



ing a similar approach. They fit a power-law t o the unmagnified HST photometry f rom the 
CASTLES survey, using the lensing models of iPeng et all (120061 ) . In comparison to IGPLIO 
we observe that our luminosity estimates are, on average, 0.20 ±0.05 dex smaller. The offset 
is likely caused by a combination of the differences in the lensing models, in the prescription 
used to deal with the flux ratio anomalies, and in the use of t he AG N SED template of 



Assef et all (120101 ) instead of the power-law fits of IPeng et all (120061 ). We note that this 
offset translates to 0.1 dex in BH mass, well below the uncertainties we estimate for our SE 
mass measurements in the next section. We also note that our conclusions are unaltered if 



we replace our 5100A continuum luminosity estimates with those of iGPLlOl for all objects 
where this is possible. 



http : //nedwww. ipac . caltech . edu/ 
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3.4. Black Hole Mass Estimates 



The width of a given broad emission line in a Type 1 AGN is primarily caused by the 
gravitational attraction of the supermassive black hole on the gas in the broad line region 
(BLR). Hence, the mass of the black hole, M B h, can be estimated from virial assumptions 
by 

Rblr(Av) 2 



Mi 



BH 



/ 



G 



(2) 



were Av is the velocity dispersion of the BLR gas, estimated from the width of the broad 
emission line, G is the gravitational constant and -Rblr is the distance from the black hole 
to the BLR. The factor / is a scale factor of o rder unity that de pends on the structure, 



kinematics and inclination of the BLR (see, e.g.. lCollin et al.ll2006l and references therein). 
The term R BLR (Av) 2 /G is usually referred to as the virial product (VP) and encapsulates 
all the observable quantities for a single objec t. The radius of the BLR can only be mea- 



sured through reverberation mapping (see, e.g.. |Peterson et al.ll2004), but has been shown to 



corre l ate well with t he continuum luminosity (see, e.g.. iKaspi et al 
20091 : IZu et aDboioh . 



2005; Bentz et al. 2006 



For the broad hydr ogen emissio n lines we estimate the BLR radius using the -Rblr — 
AL A (5100A) relation of iBentz et all (120091 ) . which was calibrated using a large sample of 
RM AGNs. The / factor of equation ([2]) depends on the characterization of the line width, 
generally either the FWHM or the line dispersion, 07, as well as on the emission line being 
used. For estimating Mbh fro m the width of the H/3 broad line, we use the / factor cali- 
brations of Collin et al. (l2006h for the FWHM and for 07. While for 07 a unique / factor of 
3.85 for all AGNs suffices, ICollin et al.l ( 120061 ) argued that / is strongly dependent on the 
line profile shape for FWHM-based estimates, where the shape was quantified as the ratio 
between the FWHM and 07. We choose, however, to use the best- fit fixed f factor of 1. 17 for 
FWHM instead of the line-shape dependent calibrations because iDenney et al.l (120091 ) have 
sho wn that 07 is affect ed by blending with other emission lines, making the correlation found 
by ICollin et al.l ( 120061 ) hard to interpret. For Ha there is no equivalent calibration of the 
/-factor, so we cannot direct l y esti mate the black hole masses. Instead, we use the relation 
determined by I Greene fc Hoi (120051 ) between the FWHM of Ha and H/3, 



FWHM H/3 = (1.07 ± 0.07) x 10 3 



FWHMh c 
10 3 km s~ 



(1.03±0.03) 



km s 



(3) 



to estimate the H/3 FWHM and then estimate Mbh (Ha) using the same /—factor and 
-Rblr — L relation as for Mbh(H/3). Unfortunately, there is no equivalent transformation for 
07, so we cannot use this measurement to estimate the mass of the black hole from Ha. 
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Combining equations © and (j3J) with the .Rblr — ALa(5100A) relation of iBentz et al. 



fl2009h we get 

M BH(Ha) . 7 .- 1o6/ (~)-(«)-m , m 

where in equation (j3J) A^h^ can be either the line dispersion or the FWHM. Because equation 
(jSJ) is fully dependent upon the scaling relations for H/3, the / factor in it is the same as for 
FWHMh/? in equation (j4]). Table |5] shows our BH mass estimates based on Ha and H/3 for 
all objects in the sample. 



For the UV/optical spectra we use the empirical Mbh calibrations of IVP06I for the Civ 
broad emission line, given by 

iu (n \ m« f_^cw V / AL A (1350A) V 53 
Mbh(Civ) = 10 — — — — — — M Q , (6) 



10 3 km s- 1 J \ 10 44 erg s" 1 J 

where Av is either FWHM or ai, and k = 6.66 ± 0.01 or 6.73 ± 0.01, respectively, for 
these line-width characterizations. The constant k implicitly contains the / factor, which is 
assumed to be a constant for all objects. Whenever possible, we use the observed 1350A flux 
to determine the continuum luminosity. Unfortunately 1350A is not within the observed 
wavelength range of all the UV/optical spectra we use. In these cases w e estimate the 
continuum luminosity at 1350A using the observed flux at 1450A, as IVP06l have shown L\ 



at these wavelengths to be equivalent. We list our C IV BH mass estimates in Table |5]for both 
prescriptions used to measure the widths of C IV. As expected, masses determined from the 
FWHM are highly consistent for both prescriptions, with a mean difference of 0.04 dex and 
a scatter of 0.02 dex, with the average prescription B based mass estimates being smaller. 
The agreement is much worse for 07, with a mean difference of 0.23 dex, in the sense that B 
is smaller, and a scatter of 0.18 dex. 

We estimate the uncertainties in our BH mass estimates by propagating the errors in 
the velocity widths and in the continuum luminosities. For masses based on the width of 
the broad Hydroge n lines, we also propagate the uncertainties in the /-factor and in -Rblr- 



Collin et al.l (120061 ) determined that the uncertainty in / when using 07 is 30%, whi le that in 



FWH M is 43%. For -Rblr we assume the intrinsic scatter of 0.11 dex estimated by [Peterson 



( 120101 ) for the radius-luminosity relation. Adding the uncertainties in / and -Rblr is not 
possible for the C IV estimates of the BH masses. Instead, we add the measurement errors 
and the in trinsic scatter between C IV and RM BH mass estimates in quadrature. Using the 



sample of IVP06I . we estimate intrinsic scatters of 0.32 dex and 0.28 dex for FWHM and ai 
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respectively. IVP06l found that the total scatter, including measurement errors, was 0.32 dex 
for both line- width characterizations of Civ, showing that the intrinsic scatter dominates 
over measurement errors, especially for FWHM estimates. 



4. Biases in C IV Black Hole Mass Estimates 

In this section we use the sample described in §[2] to study biases in the C IV black hole 
mass estimates. We first compare how the mass depends on the characterization of the C IV 
line-width, and then we proceed to compare these rest-frame UV estimates to those based 
on the Ha and H/3 emission lines. In the next section we will compare our results with those 
of other studies on the relations between C IV and H/3 BH estimated masses. 



4.1. Comparison of FWHM and o\ Derived Masses 

Given that we have measured both FWHM and o\ for Civ in all our objects, the 
simplest test we can perform is to determine if there are any biases between them as BH 
mass estimators. Both measurements have advantages, and some co ntention exists in the 



literature as to which constitutes a more reliable mass estimator (see iPeterson et al.l 12004 . 
and references therein). 

Figure compares the C iv-based BH masses determined for both line-width estimates 
and for the two continuum and line blending prescriptions A and B, respectively. A clear 
bias is observed for both prescriptions, where most objects have a lower estimated BH mass 
if we use o\ instead of the FWHM. The bias for prescription A (B) width measurements 
seems to be well represented by a constant offset of K = 0.13 ±0.06 dex (0.24 ±0.07 dex) or, 
equivalently, a factor of 1.3 (1.7). We fit for K while simultaneously fitting for the intrinsic 
scatter between the two mass estimators by adding a scatter 5" in quadrature to the error of 
each logarithmic mass difference. Note that the logarithmic mass difference does not depend 
on the continuum luminosity or the intrinsic scatter with respect to the RM estimates. In 
practice we maximize the likelihood 



C = ((a 2 ) + S 2 ) 



-1/2 



' N 

nw+ 5 T va 



where of is the variance due to measurement errors in the logarithmic mass difference of 
object i and (a 2 ) is its average over all objects. We exclude objects for which we consider the 
C iv-based BH mass estimates to be lower bounds due to absorption. The leading factor in 
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equation (JTj) is a logarithmic prior on the overall dispersion. The best fit scatter is similar for 
both prescriptions, with case A line- widths producing S = 0.16 dex while case B ones have 
S = 0.19 dex. Since the logarithmic mass difference only depends on the line- widths and 
not on the continuum luminosities, the constant BH mass offsets K can also be expressed as 
an offset between the line-width characterizations. As such, these values imply an offset of 
0.10 ± 0.03 dex (0.16 ± 0.04 dex) between the FWHM and o\ line-width characterizations of 
Civ for prescription A (B). 

It is not surprising that prescription A provides a smaller offset between BH masses 
obt ained from the FWHM and o\ of Civ, as this prescription is modeled after that used 



by IVP06l . who used their measurements to determine equation ([6]). However, given the 



similarity, the presence of a n on-zero offset for prescription A is somewhat puzzling. If we 



examine the sample of IVP06I . the scatter is larger, 0.2 dex, and there is no offset (—0.02 ± 
0.03 dex), although the lack of an offset is by definition small since both mass estimators 
were calibrated against the same RM data set. 

The large overlap in the mass and continuum luminosity ranges of our sample and that 



of IVP06I suggest that dependence on a secondary parameter is unlikely. Furthermore, we 
do not see any correlation of this bias with BH mass, continuum luminosity or Eddington 
ratio. There is also no correlation with redshift, suggesting that it is unlikely to be an 
evolutionary trend. The only other major difference between the samples is lensing by 
foreground galaxies. This, however, is very unlikely to cause such an effect, as quasars are 
quite compact and strong lensing affects the whole object. Microlensing by the stars in the 
foreground galaxy could in principle distort the shape of the C IV broad emission lines due to 
the spatial dependence of their velocity structure, but this is very unlikely for two reasons. 
First, the width of Civ is typically well below 10,000 km/s, constraining the location of 
the gas to a distance greater than > 10 3 Schwarzschild radii [Rs) from the bl ack hole, while 



micro lensing is only observed to have significant effects on scales below 100 Rs (IMorgan et al. 



20101 ). Second, the gas moving at the highest velocities is expected to be closest to the black 
hole, so microlensing would tend to magnify the wings of the line more than the core, and 
hence producing the inverse of the effect we see by making 07 too large rather than too small 
compared to the FWHM. While microlensing can also produce regions of demagnification in 
the source plane, these are of very large spatial extent, and so it is unlikely to see significant 
magnification variations across the BLR. 

It is likely then that other minor differenc es in the meth o d we use to measure 07 as com- 



pared to lVP06l give rise to the remaining bias. iDenney et all ( 120091 ) showed that estimates of 



<j\ depend on the exact prescription used for the line-width measurement and the segregation 
of blended emission for H/3. Our investigation shows that this may be the case for Civ as 
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well (see also iFine et al.ll2010l ). However, the remarkably low scatter in Figure |6] suggests 
that if a i is measured in a self-consistent manner it can be as accurate as the FWHM for 
estimating BH masses, but the calibration will depend on the exact prescription. In the 
next section we will explore the reliability of the C IV FWHM and o\ BH mass estimates by 
comparing them to those based on Ha and H/3. 



4.2. C IV compared to Ha and H/3 

Figures [7J and [S] compare the mass estimates based on the Ha and H/3 lines to those 
based on the width of Civ. We only show here (and for the rest of the figures) UV BH 
masses based on the prescription B width measurements of Civ. The FWHM based BH 
masses are almost equal for prescriptions A and B (see §3.4|) . but they show a systematic 
offset for the ai estimates (see §4.ip . We adopt the prescription B masses for the rest of this 
section, but our conclusions are unaltered if we instead use prescription A measurements. 
We have made the assumption that the Civ FWHM mass estimates are unbiased, and so 
those obtained from the prescription B o\ measurement of Civ have been shifted by the 
systematic offset of 0.24 dex derived in the previous section. 

We measure no significant offset between the Civ-based and Ha- or H/3-based masses 
when using only objects with C IV line-widths that are not lower bounds and have reliable 
Balmer line widths (group I). We find best fit offsets of -0.12±0.15, -0.11±0.16, -0.15±0.16 
and —0.19 ± 0.18 dex for panels a), b), c) and d), respectively, of Figure [HJ with residual 
scatter of 0.30, 0.23, 0.46 and 0.38 dex. Including the objects with group II Ha and H/3 
line-width estimates does not change this conclusion, with best fit offsets of —0.05 ± 0.14, 
—0.13 ± 0.13, —0.07 ± 0.15 and —0.15 ± 0.14 dex, respectively, with measured scatters of 
0.36, 0.33, 0.46 and 0.41 dex. The lack of offsets confirms our assumption that Civ FWHM 
BH masses are unbiased and that only those based on o\ need to be corrected. The constant 
offset fits yield \ 2 P er degree of freedom (xl) values of 0.6, 0.5, 1.4 and 1.1 for panels a), b), 
c) and d) of Figure [7] when using only the solid symbols. The scatter in each panel of Figure 
[7J is largely consistent with the estimated uncertainties, although the errors in the C IV ai 
masses may be slightly overestimated. We find no evidence based on the \ 2 statistic that a 
slope different from unity is required to describe the relation between the logarithms of the 
BH masses (Figure [7]) , independent of whether we include the group II Balmer line- width 
measurements. 

We next investigate if the residuals between the Civ and Balmer line masses are cor- 
related with any other observables. Figures |9] — [15] show the residuals as a function of the 
1350A and 5100A continuum luminosities, redshift, Eddington ratio, blueshift of the Civ 
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line, asymmetry of C IV (parametrized by the ratio of the widths red and blue of the cen- 
troid), and the ratio of the UV and optical continuum luminosities. Table [6] summarizes the 
significance of the correlations based on their Spearman rank-order coefficients. Only the 
correlation with the ratio of the rest-frame optical and UV continuum luminosities is signif- 
icant (Figure [T5l) . Figure [161 compares the Civ and Balmer line derived BH masses after 
rescaling the C IV masses using the best fit correlation determined from the corresponding 
panel in Figure [15j We applied corrections of the form 

logM^(Civ) = logM™(Civ) - b - a\og XLx{135 ° A) . (8) 

S AL A (5100 A) V ; 

where the coefficients a and b are listed in Table UJ For completeness, this Table also shows 
the coefficients obtained when using the prescription A line- widths of Civ, which are of 
similar magnitude and significance. Note that the uncertainties given for these coefficients 
have been determined after rescaling the errors such that the best fit has xt = !• 

The agreement between the rest-frame UV and optical BH mass estimates after applying 
this correction is remarkable, and the scatter of objects with group I and non-lower bound 
line-widths has decreased from 0.30 to 0.11, 0.23 to 0.10, 0.46 to 0.25 and 0.38 to 0.22 dex for 
panels a-d of Figures [7J and [HJ respectively. We find that the lowest scatter is between the BH 
masses estimated from the o\ of C IV <j\ and the FWHM of either Balmer line. This supports 
our conclusion in the previous section that o\ Civ BH masses have small random errors, 
even if their systematic errors may be much larger than those of the FWHM estimates due 
to blending of emission lines. Such a small scatter places strong constraints on the strength 
of a possible correlation between the mass residuals and any tertiary parameter. We find 
again that a slope different from unity is not required to describe the relation between the 
logarithm of the C IV and Balmer line BH masses. 

Since BH mass estimates generally scale as £s.v 2 L x l 2 (eqns. [4] and [6]), a naive interpre- 

1/2 1 /2 

tation of the reduced scatter is that we have simply shifted from showing L ' , vs. L ' , to 

-LooUA O-LUU-A. 

^5iooA vs ' ^5iooA' r ^^ Le ^ es ^ ^ correction is statistically different from simply replacing L 135Q ^ 
by £ 510 oA ^ l-2c, so it is not simply swapping the luminosities. More importantly, even 
if the slope was exactly a = 1/2, it reveals the crucial point that a significant 
fraction of any problems in reconciling C IV and Balmer line estimates of BH 
masses is due to the estimates of the continuum luminosities rather than any 
properties of either line. 

There are 3 potential causes for a correlation of the mass ratio with the ratio of the 
continuum luminosities: i) obscuration, ii) host contamination and iii) non-universal AGN 
SEDs. Unfortunately, our analysis does not allow us to determine which BH mass estimate is 
more accurate. Extinction will reduce the rest-frame UV continuum luminosity while having 
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little effect on the rest-frame optical luminosity. Conversely, host contamination will raise 
the optical luminosity while leaving the UV unchanged, as galaxies are typically brighter 
in the optical than in the UV. With respect to case iii), the radius of the BLR is really 
determined by the flux of the ionizing continuum (A < 912A). The -Rblr — L relations used 
to construct equations ©, © and implicitly assume a universal SED for all quasars, 
as they imply that the ionizing continuum can be uniquely predicted from the continuum 
luminosity at longer wavelengths. This approximation is likely to be better for the rest- 
frame UV continuum than for the optical. All three cases discussed would produce a slope 
of a ~ 0.5 in equation (jSJ), simply representing the luminosity power indices in equations 
©5 © and (jSj). This is generally shallower than the observed slope but within 2a of the 
best-fit relations. A larger sample is needed to fully determine if the slope of this correlation 
is statistically different from a ~ 0.5. We note that in order to create a slope larger than 
0.5, it would be necessary for the velocity widths of the quasar broad lines to be dependent 
on the ratio of the continuum luminosities. There is some evidence that the inclination 
angle of the accretion d isk with respect to th e line of sight may correlate with both the 



SED of the contin uum (Galla gher et al. 



the broad H/3 line (IWills fc Browne 



200 



1986 



and references therein) and the FWHM o f 



Wills fc Brotherton 1995 



Jarvis fc McLure 



2006), 



althoug h no such correlati on is observed for the FWHM of C IV (iVestergaard et all 120001 . 
but see iDecarli et al.1 120081 ) . Accretion disk inclination corrections, however, would act in 
the opposite sense to the observed correlation and hence cannot be responsible for a slope 
in excess of 0.5 — disks with higher inclination angles (closer to edge-on) would appear to 
have higher FWHM of Hf3 and bluer continua for a fixed "true" BH mass (i.e. not estimated 
from spectral features; see iGallagher et al.ll2005l . and references therein). 



Our sample is likely representative of observations of the general quasar population in 
terms of reddening and host contamination. It could, in principle, have a larger typical 
reddening due to additional obscuration by dust associated with the lens, but this is unlikely 
to be important for our sample. R eddening by the lens galaxy will typically vary between 
quasar images. iFalco et al.1 (119991 ) studied most of the objects in our sample and found 



that only two of them showed significant differential reddening: SBS0909+523 (AE(B — 
V) = 0.2 mag for image B with respect to A, see also Appendix A) and Q2237+0305 
(AE(B — V) = 0.18 and 0.17 mag for images C and D with respect to A). Small but non- 
zero differential reddening was also detected for three other lenses (HE1104-1805, H1413+117 
and B1422+231). The lensed quasars SDSS0246-0825, HS0810+2554, FB Q1633+3134 and 
SDSS1138+0314 were not part of the sample studied by lFalco et al.1 (Il999l ). We studied the 
latter object in §2.11 and concluded images B and C did not show evidence for differential 
reddening between them, but there is no information in this regard for the other three 
quasars. Lensing can also alter host contamination in the quasar observations as compared 
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to an unlensed case. The exact amount of host contamination depends on the size of the 
PSF and aperture used, the m orphology of the lens and the surface bri ghtness profile of the 



quasar's host galaxy (see, e.g., iKochanek et al.ll2001t iRoss et al.ll2009l ). however the zeroth 



order effect is to not alter the amount of host contamination compared to an unlensed quasar. 

While we have shown that the dominant source of scatter in the comparison between the 
BH mass estimates based on C IV and the Balmer lines is due to the continuum luminosities, 
we still wish to assess the relation between the widths of the different emission lines used. 
Figure [T7] shows the comparison between the C IV and Balmer line widths. Note that we 
do not show measurements for which we only have lower bounds on the C IV width due to 
absorption. The best agreement is between o\ of Civ and FWHM of H/3, which is expected 
given that these measurements also give the lowest scatter in the BH mass estimates, however 
a generally good agreement is also observed in all panels. We remind the reader, however, 
that the corrections we found between the BH mass estimates residuals and the ratio of the 
continuum luminosities did not have a slope of 0.5. This implies that the ratio of the line- 
widths may have a dependence on the luminosity ratio, with a power given by the excess of 
the slope from 0.5. This could be a source of additional scatter in Figure [T71 and so, instead 
of comparing the line-widths directly, we also compare them after applying a correction 
based on the continuum luminosity estimates, of the form 

, a tt \ 1 a TT \ (a -0.53), AL A (1350 A) r _ 3l AL A (5100A) 

log Av(KP or Ha) corr = log Av(W or Ha) + - log ^ + 5x 10~ 3 log J . 

(9) 

as shown in Figure [T51 Note that since we don't know the origin of the corrections, applying 
it to H/3 rather than C IV is a completely arbitrary decision made for display purposes. The 
agreement is now better and a correlation between the measurements is clear, suggesting 
that the widths of both lines are equally good tracers of BH mass. We have quantified 
the correlation between the different line- widths in Table E] using the Spearman rank-order 
coefficient, as we did for the BH mass residuals. After applying the correction from equation 
09]), we find positive correlations which are typically statistically significant between the Civ 
and Balmer line-widths. A weak anti-correlation, however, is measured between the FWHMs 
of Civ and H/3, but it is not statically significant. 



5. Comparison with Other Studies 

In §4.2l we used a sample of lensed quasars to compare BH masses based on observations 
of the C IV emission line and of the Balmer lines Ha and H/3. We found that the agreement 
between the rest-frame UV and rest-frame optical based BH masses is reasonably good. 
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We also found that this agreement is even better once we apply an empirically determined 
correction based on the ratio of the 1350A and 5100A continuum luminosities. There have 
been a number of previous studies that have explored the relative accuracy of BH masses 
based on C IV and H/3, and they have reached both similar and opposite conclusions. 



The studies of |VP06| . which we have discussed previously, and iDietrich fc Hamann 



(120041 ) found that C IV derived BH masses are consistent with those o btained from th e widt h 
of H/3, and h e nce co n stitute a valid r eplace ment as a mass estimator. IShemmer et al.l (120041 ) , 
Netzer et al.l (120071). ISulentic et al.l (120071) and IDietrich et al.l (120091 ). however, reached op- 
posite conclusions. IShemmer et al.l ( 20041 ) concluded that BH masses derived from Civ were 
poorly matched to those obtained from H/3 and could be systematically different. They 
showed that for a sample of narrow-line Seyfert 1 galaxies, C IV based BH masses are larger 
by an average factor of ~ 3 with respect to those obtained from H/3. IDietrich et al.l (120091 ) 
also found a large disagreement between the two estimates of the BH mass, but they found 
that using Civ tends to underestimate the BH masses by a factor of ~ 1.7, although the 
significance of t his re sult is limited by the small number of objects (9) in their sample. While 
Sulentic et al.l (120071 ) also found significant disagreement, they argue t hat the magn i tude o f 
the offset depends on the spectroscopic characteristics of the quasar. iNetzer et al.l ( 120071 ). 
on the other hand, found no significant offset between the mass estimates, but argued there 
was also no discernible correlation between them. It is likely that many of the differences 
between the results of these studies are due to the use of different -Rblr ~ L calibrations, 
different /-factors, different prescriptions for measuring line widths, limited mass ranges and 
data quality. 

Here we take their measurements, where possible, and make estimates of the BH masses 
using equations (@|, (J5]) and (Q. We caution the reader, however, that we are not redoing 
the line-width and continuum luminosity measurements in a consistent manner and that 
this may b e a significant sourc e of additional scatter. W e used all 21, 15 and 9 sources 
from IvPpd . INetzer et all toOl\ ) and bietrich et all koO$ ) with H/3 and C IV line FWHM 
and continuum luminosity measurements. We could not use 10, 29 and 1 sources from these 
studie s that lack either or both line widths, or any of the sources from IDietrich fc Hamann 
( 2004 ). which lac k measurement s of th e 5 100 A continuum luminosity. We also could not 
use the sources of ISulentic et al.l ( 120071 ). as they only report narrow-component subtracted 
C IV widths, which are not compatible with the rest of the measuremen t s we discuss. We 
note that the 29 sources we could not u se from the study of INetzer et al.l (120071 ) also belong 
to the sample of IShemmer et al.l (120041 ). for which the Civ line widths and UV continuum 
luminosities are not reported. 



The left panel of Figure [19] compares the C IV and H/3 BH masses derived for all these 
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objects along with those in our sample. A clear correlation is observed for the complete 
ensemble of objects, albeit with a considerable scatter of 0.41 dex. The scatter is comparable 
with the 0.46 dex we find for our sample of C IV and H/3 FWHM-based BH masses (see §4.2p . 
A Spearman rank-order coefficient analysis returns r s = 0.79 with a probability that both 
mass estimates are uncorrelated of P ran = 2 x 10~ 12 . A linear fit to the left panel of Figure 
[T9l of the form 

Mbh(Civ) M bh (H/3) 
log 3 = m log ^ + n (10) 

& 10 8 M Q & 10 8 M K J 

returns a best-fit slope of m = 0.89±0.08 and intercept of n — — 0.09±0.08 (the measurement 
uncertainties were scaled to make xl = 1 before determining the uncertainties in the fit 
parameters). If we plot the residuals between the two BH mass estimates we find, just as 
in §4.21 that a significant correlation is observed with the ratio of the UV and continuum 
luminosities (Figure [20]) . but not with BH mass, redshift, Eddington ratio or the continuum 
luminosity at 5100A (all shown in Figure [2T|) . or with the continuum luminosity at 1350A 
(not shown). The best-fit linear relation to the correlation between BH mass residuals and 
the ratio of the continuum luminosities, shown in Tabled has a slope of a = 0.82 ± 0.18 and 
an intercept of b = —0.40 ± 0.07. While the slope is consistent with the value obtained for 
our sample alone (a = 0.86 ± 0.25, b = —0.23 ± 0.12), the intercept differs by approximately 
0.2 dex (approximately 1.5a). The offset is likely produced by the different prescriptions 
used to measure the width of the broad emission lines. Figure [19] also compares the C IV 
and H/3 derived masses after correcting for this correlation by applying equation (jSJ) (see 14.21 
for details). While the strength of the correlation has not increased substantially (r 8 = 0.80, 
-Pran = 6.4 x 10 -13 ), the scatter has decreased from 0.41 dex to 0.34 dex. This change is 
significantly more modest than that found for our sample of lensed quasars, but this is likely 
due to the inhomogeneous prescriptions used to measure the width of the emission lines. A 
linear fit of the form of equation (TTUT) to the relation between the BH mass estimates after 
applying the correction returns very similar parameters as before, with a best-fit slope of 
m = 0.88 ± 0.07 and intercept of n — 0.06 ± 0.07. Note that the measurement errors have 
again been scaled to make \ 2 V = \ before estimating the uncertainties in the fit parameters. 
Given the inhomogeneity of the measurements used, the relatively small number statistics 
of the sample, their typically large error bars and the likely intrinsic dispersion, however, we 
cannot currently determine whether the deviation from a slope of unity is significant or not. 

As argued before, the inhomogeneity of the measurements can be a very significant 
source of scatter in the comparisons discussed above. We have shown in the previous section 
that with homogeneously analyzed, high S/N spectra, the difference between the line widths 
is not the dominant source of scatter in the comparison between C IV and Balmer line-based 
BH mass estimates. We would still like to assess whether the Civ and H/3 line widths are 
correlated in this combined sample. Figure [22] shows this comparison for all the objects 
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used in this section with and without applying the correction based on equation (Q. While 
the scatter in Figure [22] is large, there is still a statistically significant (99%) correlation 
between the measu red line- widths (see Table [S]). Most of the scatter is due to the sample of 
Netzer et al.l ( 120071 ). Upon inspection of the SDSS spectra used for that study, we find that 
almost all the outliers correspond to low S/N spectra. Given that the Civ line is typically 
very complex, this can be a major source of uncertainty. 

As an experiment, we obtained higher S/N spectra for one of the outliers, SDSS1151+0340. 
It has the third most discrepant line- width ratio in the sense that the C IV line is too narrow 
compared to the Ba lmer lines. We obtained two independent spectra, one with OSMOS 



(Martini et al 



( lOke fc Gunn 



20111) at the MDM 2.4m telescope and one with the Double Spectrograph 



1982) at the Palomar 200-inch telescope. Due to poor weather conditions and 



aperture size, only the Double Spectrograph observations yielded a higher S/N spectrum 
than that of SDSS. All three spectra of SDSS1151+0340 are shown in Figure [231 The spec- 
trum obtained with Double Spectrograph reveals that there is significant absorption near the 
Civ line, with two clear absorption troughs. These can be seen in the lower S/N spectra, 
but are difficult to distinguish from the noise. Due to the very substantial absorption, it is 
not possible t o reli ably measure the width of the Civ line, and the width measurement of 
Netzer et al.l ( 120071 ) should only be considered as a lower bound. While the SDSS spectrum 
of this source has the lowest continuum S/N in their sample (S/N = 1.5), it is comparable 
to many of their other sources. The average continuum S/N of the SDSS spectra is only 6.7, 
with all objects having a lower S/N than any optical spectra in our lensed quasar sample. 
In particular, the second largest outlier in their sample also has the second lowest S/N of 
4.1. 



While our example comes from iNetzer et al.l (120071 ) , low S/ N spectra are also present 
in all the additional samples we consider. If we eliminate objects with continuum S/N < 10 
in the vicinity of C IV, the statistical correlation between H/3 and C IV line widths increases 
dramatically. The bottom panels of Figure [22] show the comparison of the line- widths in the 
absence of these objects, and a clear correlation is observed between the C IV and H/3 FWHM 
measurements, regardless of whether we apply the continuum luminosity based correction 
of equation l^fy. These correlations are more statistically significant by about two orders of 
magnitude than when including the low S/N spectra, with a probability of not being real of 
~ 5 x 10~ 4 (see Table [8]for details). This suggests that the width of C IV is as good a tracer 
of BH mass as the width of H/3, with the caveat that high S/N spectra of the rest-frame UV 
region are fundamental to accurately model the structure of the C IV emission line. 
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Conclusions 



Using a sample of high-redshift gravitationally lensed quasars observed spectroscopically 
in the UV/optical and NIR, we have studied the agreement between single-epoch BH mass 
estimators based on the ClV, H/3 and Ha broad emission lines. Our sample consists of 12 
lensed quasars observed with HST by the CAST LES pro ject. In particular, we have used 
the sample of NIR spectroscopic observations by iGPLlOl as a starting point and improved 
on it by (i) adding new NIR observations for 3 objects (SDSS1138+0314, SBS0909+253 and 
HS0810+2554), (ii) adding high S/N, uniformly analyzed, optical spectroscopic observations 
for all targets, and (iii) adding the missing rest-frame \L\(5100A) luminosity estimates for 
SDSS0246-0852, HS0810+2554 and Q2237+030. 

We described in detail all the methods we used to measure velocity widths and their 
uncertainties, rest-frame continuum luminosities and to estimate the BH mass of each quasar 
using the H/3, Ha and Civ emission lines. We first compared the Civ BH mas s estim ates 
based on the FWHM and o\ line-width characterizations and the calibration of IVP06I and 
found that, for our sample, the o\ based BH masses are systematically underestimated with 
respect to the FWHM-based ones by 0.13 ± 0.06 dex if using prescri ption A and 0.24 ± 
0.07 dex if using prescription B. A similar offset is not observed in the IVP06l data set. The 
difference probably arises from our treatment of the blending of the broad Civ emission 
line with the nearby broad Hell A1640 and Fell emission redward of C IV, which is partly 
confirmed by the lo wer difference f o und f or the prescription A measurements. This adds 
to the arguments in iDenney et al.l (120091 ) that <j\ is not universally reliable for SE mass 
estimates in the presence of blending, as the results obtained are highly dependent on the 
exact prescription used for the line characterization. When comparing with BH masses 
derived from the Ha and H/3 broad emission lines, we find that Civ FWHM based BH 
masses are not biased, reinforcing the conclusion that the bias is in the <j\ estimates. We 
note, however, that the scatter between Civ FWHM and o\ derived masses is relatively 
small, suggesting that if a consistent prescription for measuring o\ is applied, o\ would be 
at least as accurate as FWHM. This is important because o\ measurements are significantly 
more relia ble for complex line profile shapes and i n the presence of narrow-line component 
residuals ( [Peterson et al.l 12004 ; IDenney et al.ll2009[ ). 



We then compared the C IV and Balmer line BH mass estimates. After offsetting the C IV 
<7; masses to agree with the FWHM estimates, we find there is no significant offset between 
Civ and either Balmer line BH mass estimates. Averaged over the 4 possible Civ/Balmer 
line mass comparisons (see, for example, Figure [7]), the offset is —0.15 ± 0.17 dex and the 
scatter is 0.35 dex. Note that the error in the mean offset corresponds to the average of 
the errors of the four estimates, which is representative given that the estimates are not 
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truly independen t. The scatter of 0.35 dex is very close to the scatter of 0.34 dex found by 
Shen et al.l (120081 ) between Mg II and C IV FWHM based BH mass estimates, and significantly 
larger than the scatter of 0.22 dex they found between Mgll and H/3 FWHM based BH 
masses. 

We find that the residuals between the C IV and H/3 and Ho; based mass estimates are 
not strongly correlated with the UV or optical continuum luminosities, redshift or Eddington 
ratio, but we find a strong dependence on the ratio of the UV to optical continuum lumi- 
nosities. If we correct for this color dependence, the agreement between the C IV and Balmer 
line estimates is remarkably good, with an average scatter of 0.18 dex, almost a factor of 2 
smaller. We find the scatter is smallest — approximately 0.1 dex — when using the H/3 line 
and the <j[ characterization of C IV rather than its FWHM. This observed correlation could 
be caused by i) reddening, ii) host contamination, or iii) an object-dependent SED shape. 
The slope we observe is somewhat steeper than that expected in any of these cases, and may 
suggest a luminosity component to the line-width characterization of the broad emission 
lines. A larger sample is needed to accurately determine the slope of this correlation and 
determine its nature with certainty. More generally, the comparison shows that many of 
the problems in comparing C IV and Balmer line BH mass estimates are associated with the 
continuum luminosities rather than any potential physical complexities with the C IV lines. 
When we compare the line-widths directly instead of the BH masses, we find that the width 
of Civ is well correlated with those of the Balmer lines once the correction based on the 
ratio of the continuum luminosities is applied. 

Our conclusions are un change d if we add 45 addit ional , but heteroKeneously analyzed, 
C IV and H/3 estimates from lvPod . iNetzer et all J2007f ) and bietrich et al.l (£3). We used 
the published FWHM of both emission lines and rest-frame UV and optical continuum 
luminosities of these sources, but the mass calibrations used for our sample. There is a clear 
linear correlation between the BH mass estimates, and the residuals are again correlated 
with the ratio of the continuum luminosities. The residuals are not correlated with either 
continuum luminosity alone, redshift, BH mass or Eddington ratio. We also find for this 
heterogeneous sample that the width of Civ is well correlated with that of H/3, particularly 
after we eliminate the objects with low S/N Civ spectra. Relatively high S/N spectra are 
essential to obtaining accurate line widths. 

In summary, our results show that C IV is a good BH mass estimator but with small 
prescription-dependent offsets. The correlation of the mass residuals with the continuum 
slope could be a bias in either or both of the estimators. Determining the "blame" would 
require an independent mass estim ate, but its exist ence should not be a surprise given that 
quasar SEDs are not universal (e.g. jYip et al.ll2004j ). More generally, unless we are to believe 



-28 - 



that all properties of AGN are determined by a single quantity, the black hole mass, both 
single-epoch mass estimates and reverberation-mapping radius estimates must depend on 
additional parameters. That the black hole mass seems to dominate is convenient, but the 
excess scatter in mass and radius estimates beyond the measurement uncertainties requires 
either that the error estimates are incorrect or is evidence for additional para meters. One 
possibility is that radiation pressure plays a significant role (IMarconi et al.ll2008l ) and it could 
easily affect different lines in different ways. While there has been consid erable recent effort 
to expand the range of black hole m asses included in these studies (e.g., iKaspi et al.l 120071 ; 



Bentz et al.ll2009t iBotti et al.l |2010| ). it is equally important to expand the range in other 



physical parameters such as spectral shape and Eddington ratio in order to better search for 
these additional correlations. 
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A. Notes on Individual Objects 

In this section we discuss some details of our line-width and continuum measurements 
of individual objects. All UV/optical spectra, as well as the continuum and line-profile fits, 
are shown in Figure H LUCIFER spectra of SDSS1138+0314 and HS0810+2554 are shown 
in Figures [T] and [2J while the LIRIS spectra of SBS0909+532 are shown in Figure EJ 



HS0810+2554 — The CIVA1549 profile of HS0810 shows a small amount of absorption 
near the peak of the line. We interpolate over this region before making the line-width 
measurements and fitting the GH polynomial. Our results are consistent with or with- 
out the interpolation, as the absorption is weak and only seen near the very peak of the 
line. To fit the continuum and emission-line features that blended with the H/3 emission 
of HS0810+2554, a power-law continuum and Fell broad-emission line template were fit to 
the spectrum ba sed on the continuum re gions listed in Table H] and the rest-frame optical 



Dietrich et al. 20021 



2005 



Fe H template of iBoroson fc Green! (119921 ) from observations of I Zwl (see I Wills et al. 



1985 



for more details). Narrow [O ill] AA4959,5007 emission was then 
removed by creating a template from a two-component Gaussian fit to the [O ill] A5007 nar- 
row line and then scaling it to [O ill] A4959 based on standard emission line ratios. We could 
not remove narrow H/3 emission because such a component was not obvious in the spec- 
trurrfl After subtracting these components, the remaining broad H/3 emission was fit with 
a Gauss-Hermite polynomial, and the FWHM and line dispersion were measured from this 
fit as described in §3.11 The deblended spectrum of HS0810+2554, showing each component 
including th e GH fit, is shown in F igure |2j Our H/3 FWHM measured from the LUCIFER 
spectrum of iMogren et al.l (jin prep.l ) is consistent with that of IGPLIO . 



SBS0909+532 — We use the combined UV-optical-NIR spectrum of iMediavilla et al.l (120 lOl ) 
of images A and B of this object, based on a combination of HSTSTIS and WHT INTEGRAL 
and LIRIS observations. The UV section of the spectrum is shown in Figure H] while the NIR 
section is shown in Figure [3j The C IV profile of SBS0909+532 showed a small absorption 
trough near observed frame 3600A and we interpolated over this region before measuring 
widths. The SBS0909+532 C IV profile shape is 'peaky' with broader wings at the base, and 
our GH fitting procedure was unable to satisfactorily fit this line profile, so we estimate errors 
based on the original spectra instead of a GH polynomial fit. For this object we only measure 
the UV continuum luminosity on image A , as im age B shows clea r differe ntial reddening with 
respect to A. We note that the lPeng et al.l (120061 ) mass quoted by lGPLlOl is based on Mgll, so 
using the C IV line- width measurements given in Table El provides the first estimate of a C IV- 



GPL10 are similarly unable to isolate a narrow component in their observations of HS0810+2554 
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based black hole mass for this object. For the IR spectra of iMediavilla et al.l ( 120101 ). shown 
in Figure |3l we removed narrow-line components from the IR spectra using the [O ill] A5007 
line as a template and scaling it to the other narrow lines using standard emission-line ratios 
between lines of the same atomic species and basing the strength of the Balmer narrow lines 
on the ratio of [Om] A5007/H/3 determined by inspection. We are not as confident in our 
narrow-line subtraction for this object as for the others because (1) we see residuals near 
the peak of H/3, and (2) the exact strength of Ha is uncertain because narrow- line emission 
remains present after subtraction. The exact level of the residuals for Ha is unclear, since 
increasing the fraction of emission by as much as a factor of 2 does not result in an obviously 
improved subtraction. In the case of H/3, the residuals are not larger than expected based 
on the S/N of the images, but for Ha we report uncertainties determined from difference 
between the widths determined with or without the narrow-line subtraction. This results 
in an Ha FWHM uncertainty several times larger than would be estimated by our Monte 
Carlo simulations. Comparable o\ uncertainties are measured using both methods, because 
the line dispersion is far less dependent on the presence of a narrow-line component (see 



Penney et al.ll2009l ). We measure the H/3 line- widths from the GH fits to the profile and the 
Ha line-widths directly from the data because the GH polynomials did not accurately fit the 
line profile. Image B may have a residual sky line peak just blueward of the H/3 narrow-line 
component. The presence of this emission has little effect on our fits, however, since we 
measure consistent line widths if we interpolate und er this e mission to remove it. Our Ha 



and H/3 widths are consistent with those reported by IGPLIO 



Q095 7+561 — We use the HSTSTIS UV spectrum of both images obtained by lGoicoechea et al. 



( 120051 ). The rather strange CivA1549 line profiles in this object may indicate that there is 
absorption and/or that the profile shapes in individual images are affected by microlensing 
from the lens galaxy. However, since there was no definite source of uncertainty to correct 
for, we simply measured the observed line widths from each spectrum. 



HE11 04-1805 — We use the EFOSC1 ESQ 3.6m telescope UV/optical spectrum of lWisotzki et al 



( 119951 ). The C IV A1549 profile shows a small amount of absorption near the peak of the line, 
similar to that of the HS0810+2554 CivA1549 profile. We therefore apply the same treat- 
ment to this line as to the HS0810+2554 profile, and find similarly consistent results with 
or without interpolation. 



PG1115+080 - Due to the severe absorption, both narrow and broad, in the CivA1549 
line profile in this object, we could not measure the CivA1549 line-width directly from the 
data. However, by masking out the absorption regions, we made a reasonable GH fit to this 
line profile, from which we measured the line widths given in Table EJ 
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SDSS1138+0314 — To measure the width of Civ and the UV continuum luminos ity we 
use the FORS1 VLT spectra of images B and C obtained by lEigenbrod et al.l (120061 ) . The 
CivA1549 line profile not only shows absorption in the blue side of the line, but is also 
particularly narrow and 'peaky' with a broad base, similar to SBS0909+523. We were unable 
to reasonably approximate the profile shape with a sixth-order Gauss-Hermite polynomial. 
However, since the S/N of this spectrum was very high (see Table [2]), we interpolated 
over the absorption with a 2nd order polynomial, measured the line width directly from 
the interpolated data, and used this interpolated spectrum and the error spectrum formed 
with equation ([1]) to derive uncertainties in the CivA1549 width measurement. Because of 
the combined effects of absorption and the narrow line profile (i.e., where the absorption 
could be masking the true width), we treat our Civ widths as lower limits. At rest-frame 
optical wavelengths, the difficulty in removing the blended narrow-line components of Ha 
and [N n] AA6548, 6583 combined with our attempt to accurately fit the emission-line peak 
(often underestimated with line profile fits) led to an overestimate of the flux between the 
Ha and NllA6583 narrow lines. This overesti mate do es not significantly affect our width 
measurements. This object was not part of the iGPLlOl sample. 



H1413+117 - This object is a BAL QSO and therefore a large portion of the CivA1549 
line profile is completely absorbed on the blue side. Hence, we adopt the C IV A 1549 width 
measured from only the red side of the line, and we consider this to be a lower limit on the 
width. 



B1422+231 - We use the LRIS Keck II UV/optical spectrum of lTonrv! (119981 ). We inter- 
polated over the two small absorption troughs near ~6875A and ~7020A before measuring 
the CivA1549 widths directly from the dat a. Our treatment of these regions did not affect 
the resulting GH fit to the data. From the IGPLlOl data, we cannot assess the reliability of 
their fit to the H/3 profile, because they plotted the H/3 spectrum of HE1104-1805 in place 
of the spectrum of B1422+231. In order to be conservative , we the refore flag B1422+231 as 



one of the objects with possible problems in the sample of IGPLIO 



FBQ1633+3134 - There is evidence for absorption in the blue side of Civ, however, it 
is not clear that a reliable interpolation could be made across this possible absorption. We 
measure the line width as is, and treat this measurement as a lower limit. 



0,2 231+030 - We use th e FORS1 VLT UV/optical spectra of images C and D obtained 



by lEigenbrod et al.l (120081 ) . We follow the same prescription as for SDSS1138+0314 and 
Q0957+561 and use an average of the Civ line-widths of each image to estimate Mbh- 
This object shows C IV A1549 absorption in the red side of the line. We interpolate over 
this absorption with a third order polynomial before measuring the line width and fitting 
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the GH polynomial to the data. The interpolation creates a peak slightly higher than that 
observed in the original spectrum, but makes for a much more symmetric line profile, which is 
more typical of the core of CivA1549 line profiles, than a linear or quadratic interpolation. 
This increase in the assumed line peak decreases our line-width measurements, but not 
significantly. 
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Table 1. Lens Magnifications and Continuum Luminosities 



Object 


Ref. 




z 


Imaj 


ie Ma 


gnification 




log 


\L X / erg 


s- 1 












A 


B 


C 


D 


(mag) 


1350A 


1450A 


5100A 


Q0142-100 


a 


2 


72 


3.3 


0.4 






16.56 


46.83 


46.76 


46 


27 


SDSS0246-0825 


a 


1 


69 


26.9 


8.9 






20.39 




44.53 


44 


59 


HS0810+2554 




1 


51 


47.2 


51.1 


13.5 


7.7 


18.72 




44.44 


44 


84 


SBS0909+532 


b 


1 


38 


1.7 


1.5 






15.18 


46.08 


46.05 


46 


31 


Q0957+561 


c 


1 


41 


3.1 


1.7 






16.51 


46.31 


46.25 


45 


79 


HE1104-1805 


d 


2 


32 


16.2 


2.3 






18.52 


46.15 


46.09 


45 


38 


PG1115+080 


a 


1 


72 


19.6 


18.7 


3.2 


4.9 


19.13 




45.47 


44 


93 


SDSS1138+0314 


e 


2 


44 


7.3 


3.7 


5.2 


6.9 


20.65 


44.83 


44.77 


44 


81 


H1413+117 


a 


2 


55 


8.2 


6.8 


6.8 


3.4 


18.05 


45.73 


45.78 


45 


63 


B1422+231 


f 


3 


62 


6.6 


8.2 


4.3 


0.3 


16.55 


46.83 


46.74 


46 


42 


FBQ1633+3134 




1 


52 


2.7 


0.7 






16.85 


45.65 


45.64 


45 


72 


Q2237+030 


9 


1 


69 


4.9 


4.3 


2.2 


4.1 


16.83 




45.53 


45 


98 



Note. - Literature UV /optical spectra o btain ed f rom the following r eferen ces: a) SPSS DR7 



Abazaiian et al. 


2009 


), b) 


(1995 


),e) 


Eigenbrod et al. 



(2010), c) 

mm,g) 



Goicoechea et a 



Eigenbrod et al. 



(120051) . d) IWisotzki et al. 



( 120081 ) . 



^Unmagnified quasar if-band magnitudes. 



Table 2. C IV Emission Line and Continuum Region Boundaries 



Prescription A Prescription B 

Blue Cont. Red Cont. Broad Line Blue Cont. Red Cont. Broad Line Res 

Object s/N* (A) (A) (A) (A) (A) (A) (A) 



Q0142-100 

V y V 7 -L -L _L V 7 V 7 


45 


5401- 


-5468 


6250- 


-6343 


5570- 


-5954 


5440- 


-5500 


5900- 


-5960 


5500- 


-5900 


2.6 


SDSS0246-0825 


17 


3980- 


-4030 


4519- 


-4586 


4028- 


-4306 


3980- 


-4030 


4250- 


-4300 


4050- 


-4250 


2.6 


HS0810+2554 


19 


3610- 


-3690 


4217- 


-4280 


3758- 


-4018 


3620- 


-3685 


3985- 


-4005 


3715- 


-3975 


15.2 


SBS0909+532-A 


71 


3420- 


-3470 


4070- 


-4120 


3564- 


-3810 


3420- 


-3470 


3820- 


-3850 


3560- 


-3820 


3.5 


SBS0909+532-B 


16 


3420- 


-3470 


4070- 


-4120 


3564- 


-3810 


3420- 


-3470 


3820- 


-3850 


3560- 


-3820 


3.5 


Q0957+561-A 


66 


3434- 


-3543 


4049- 


-4109 


3609- 


-3858 


3550- 


-3600 


3810- 


-3830 


3655- 


-3805 


3.5 


Q0957+561-B 


28 


3434- 


-3543 


4049- 


-4109 


3609- 


-3858 


3550- 


-3600 


3810- 


-3830 


3655- 


-3805 


3.5 


HE1104-1805 


221 


4731- 


-4880 


5578- 


-5661 


4971- 


-5314 


4870- 


-4910 


5300- 


-5340 


4970- 


-5290 




PG1115+080 


86 


3897- 


-4020 


4595- 


-4663 


4095- 


-4378 


4000- 


-4031 


4415- 


-4455 


4085- 


-4415 


2.2 


SDSS1138+0314-B 


78 


4904- 


-5059 


5782- 


-5868 


5153- 


-5509 


4974- 


-5020 


5520- 


-5540 


5020- 


-5520 




SDSS1138+0314-C 


47 


4904- 


-5059 


5782- 


-5868 


5153- 


-5509 


4974- 


-5020 


5520- 


-5540 


5020- 


-5520 




H1413+117 


43 


5073- 


-5233 


5981- 


-6070 


5331- 


-5698 


5140- 


-5190 


5630- 


-5680 


5210- 


-5615 


2.3 


B1422+231 


270 


6598- 


-6806 


7778- 


-7894 


6933- 


-7411 


6690- 


-6740 


7400- 


-7450 


6830- 


-7400 




FBQ1633+3134 


52 


3350- 


-3412 


4234- 


-4297 


3773- 


-4034 


3350- 


-3412 


4020- 


-4060 


3795- 


-3980 


13.0 


Q2237+030-C 


36 


3840- 


-3962 


4528- 


-4595 


4035- 


-4314 


3870- 


-3910 


4275- 


-4315 


4040- 


-4265 




Q2237+030-D 


54 


3840- 


-3962 


4528- 


-4595 


4035- 


-4314 


3870- 


-3910 


4275- 


-4315 


4040- 


-4265 





a The S/N quoted is the S/N per pixel averaged over all continuum regions listed for each object. 



Table 3. Velocity Widths 









Broad Emission Line Velocity Widths / 10 3 km s 








Object 


FWHMc iv-p-A 


FWHMc iv-P-B 


°"i,civ-p-A 


°"!,Civ-p-B 


FWHMh/s 1 " 




FWHMhq 1 " 


ff I,Ha 


Q0142-100 


5.20 ± 0.18 


4.75 ± 0.22 


3.67 ± 0.04 


3.67 ± 0.07 


2.70 ± 0.60* 




3.80 ± 0.30* 




SDSS0246-0825 


4.43 ± 0.24 


4.40 ± 0.22 


3.09 ± 0.06 


2.51 ± 0.04 


2.50 ± 0.60 




2.50 ± 0.20 




HS0810+2554 


3.68 ± 0.14 


3.53 ± 0.16 


3.40 ± 0.05 


2.85 ± 0.08 


4.40 ± 0.06 


2.21 ± 0.04 


3.80 ± 0.00 




SBS0909+532 


2.38 ± 0.08 


2.36 ± 0.08 


2.90 ± 0.04 


2.82 ± 0.05 


3.95 ± 0.17 


2.20 ± 0.05 


3.06 ± 0.34 


5.24 ± 0.04 


Q0957+561 


3.68 ± 0.25 


3.47 ± 0.08 


3.27 ± 0.09 


2.272 ± 0.007 


3.30 ± 0.90 




3.00 ± 0.20 




HE1104-1805 


6.08 ± 0.35 


5.75 ± 0.05 


3.84 ± 0.10 


2.897 ± 0.004 


3.80 ± 0.90 




4.70 ± 0.20* 




PG1115+080 


4.98 ± 0.18 


4.67 ± 0.13 


3.68 ± 0.04 


3.40 ± 0.04 


4.40 ± 0.20 




4.00 ± 0.10 




SDSS1138+0314 


>2.02 ± 0.15 


>1.99 ± 0.18 


>3.12 ± 0.04 


>2.40 ± 0.11 


3.93 ± 0.30 


2.08 ± 0.18 


2.57 ± 0.04 


1.90 ± 0.05 


H1413+117 


>2.62 ± 0.95 


>2.54 ± 0.37 


>3.78 ± 0.15 


>1.82 ± 0.07 


6.70 ± 1.90 




5.30 ± 0.80 




B1422+231 


5.81 ± 0.16 


5.56 ± 0.02 


3.70 ± 0.03 


3.321 ± 0.006 


6.10 ± 2.20* 








FBQ1633+3134 


>4.71 ± 0.18 


>4.40 ± 0.16 


>3.83 ± 0.06 


>2.20 ± 0.06 


4.60 ± 0.90* 




4.10 ± 0.70 




Q2237+030 


3.96 ± 0.18 


3.78 ± 0.12 


3.49 ± 0.07 


2.51 ± 0.05 


3.80 ± 1.40 




4.80 ± 0.60* 





tAll Ha and H^ line width measurements correspond to those in Table 1 of lGPLld. except for SDSS1138+0314, SBS0909+523 and the H/3 widths of 
HS0810+2554. 



*Group II line-widths. See ^13. II for details. 
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Table 4. NIR Emission Line and Continuum Region Boundaries 





Emission 




Blue 


Cont. 


Red Cont. 


Broad Line 


Res 


Object 


Line 


S/N a 


(A) 


(A) 


(A) 


(A) 


HS0810+2554 a 


H/3 


36 


4680- 


-4710 


5080- 


-5120 


4710- 


-4960 


8.0 


SBS0909+532-A 


H/3 


68 


11205- 


-11310 


11840- 


-11855 


11380- 


-11840 




SBS0909+532-B 


H/3 


22 


11205- 


-11310 


11840- 


-11855 


11380- 


-11840 




SBS0909+532-A 


Ha 


25 


14210- 


-14440 


17500- 


-17700 


14865- 


-16650 




SBS0909+532-B 


Ha 


18 


14210- 


-14440 


17500- 


-17700 


14865- 


-16650 




SDSS1138+0314 


H/3 


12 


16150- 


-16300 


17485- 


-17623 


16300- 


-17100 


8.0 


SDSS1138+0314 


Ha 


8 


21780- 


-21850 


23270- 


-23310 


22135- 


-23040 


8.0 



a Rest frame wavelengths are used here because line boundaries were chosen after the 
deblending procedure that transfers the spectrum into the rest frame. 



Table 5. H/3, Ha and CIV BH Mass Estimates 



Object 








log M-qyl/Mq 








FWHMciv-p-A 


FWHMciv-p-B 


(7 1 r 1 TV - 


-p.A 


07,civ-p-B 


FWHMhs 


(7 1 Y{[3 


FWHMhq 


Q0142-100 


9.59 ± 0.32 


9.51 ± 0.33 


9.36 ± 


0.29 


9.36 ± 0.29 


8.94 ± 0.30* 




9.33 ± 0.23* 


SDSS0246-0825 


8.24 ± 0.33 


8.23 ± 0.33 


7.99 ± 


0.29 


7.81 ± 0.29 


8.00 ± 0.31 




8.08 ± 0.23 


HS0810+2554 


8.03 ± 0.32 


7.99 ± 0.33 


8.02 ± 


0.29 


7 87 ± 29 


8.62 ± 0.22 


8.54 ± 0.17 


8.65 ± 0.22 


SBS0909+532 


8.51 ± 0.32 


8.51 ± 0.32 


8.76 ± 


0.29 


8.73 ± 0.29 


9.29 ± 0.23 


9.29 ± 0.17 


9.15 ± 0.24 


Q0957+561 


9.02 ± 0.33 


8.97 ± 0.32 


8.98 ± 


0.29 


8.67 ± 0.29 


8.86 ± 0.33 




8.87 ± 0.23 


HE1104-1805 


9.37 ± 0.33 


9.32 ± 0.32 


9.04 ± 


0.29 


8.79 ± 0.29 


8.77 ± 0.30 




9.05 ± 0.23* 


PG1115+080 


8.83 ± 0.32 


8.78 ± 0.32 


8.64 ± 


0.29 


8.57 ± 0.29 


8.66 ± 0.23 




8.68 ± 0.22 


SDSS1138+0314 


7.71 ± 0.33* 


7.69 ± 0.33* 


8.15 ± 


0.29* 


7.93 ± 0.29* 


8.50 ± 0.23 


8.47 ± 0.19 


8.22 ± 0.22 


H1413+117 


8.41 ± 0.45* 


8.39 ± 0.35* 


8.80 ± 


0.29* 


8.17 ± 0.29* 


9.39 ± 0.33 




9.29 ± 0.26 


B1422+231 


9.69 ± 0.32 


9.65 ± 0.32 


9.37 ± 


0.29 


9.27 ± 0.29 


9.72 ± 0.38* 






FBQ1633+3134 


8.88 ± 0.32* 


8.82 ± 0.32* 


8.77 ± 


0.29* 


8.29 ± 0.29* 


9.11 ± 0.28* 




9.11 ± 0.27 


Q2237+030 


8.67 ± 0.33 


8.63 ± 0.32 


8.63 ± 


0.29 


8.34 ± 0.29 


9.08 ± 0.39 




9.38 ± 0.25* 



Note. — All BH masses correspond to those obtained from eqns. (j4|), ([5]) and ([6]). None of the corrections discussed in § §4.11 
and 14.21 have been applied. 

*Based on group II line-width. See ^3.11 for details. 

*Should be considered as lower bound. See Appendix [A] for details. 
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Table 6. Correlations of BH Mass Residuals 



Variable 


Civ 


Balmer 


Group I Estimates 


Group I & II Estimates 




Width 


Line 


r a 


N 


^"ran 


r a 


N 


^ran 


Mbh Balmer lines 




H/3 


-0.214 


7 


0.644 


-0.183 


9 


0.637 






Ha 


0.100 


5 


0.873 


-0.048 


8 


0.911 




FWHM 


H/3 


-0.500 


7 


0.253 


-0.333 


9 


0.381 






Ha 


-0.400 


5 


0.505 


-0.143 


8 


0.736 


AL A (1350A) 




H/3 


0.357 


7 


0.432 


0.467 


9 


0.205 






Ha 


0.500 


5 


0.391 


0.571 


8 


0.139 




FWHM 


H/3 


0.143 


7 


0.760 


0.367 


9 


0.332 






Ha 


0.100 


5 


0.873 


0.429 


8 


0.289 


AL A (5100A) 




H/3 


-0.214 


7 


0.644 


-0.067 


9 


0.865 






Ha 


-0.300 


5 


0.624 


-0.048 


8 


0.911 




FWHM 


H/3 


-0.500 


7 


0.253 


-0.233 


9 


0.546 






Ha 


-0.700 


5 


0.188 


-0.214 


8 


0.610 


Rcdshift 


<ri 


H/3 


0.607 


7 


0.148 


0.517 


9 


0.154 






Ha 


0.500 


5 


0.391 


0.500 


8 


0.207 




FWHM 


H/3 


0.750 


7 


0.052 


0.583 


9 


0.099 






Ha 


0.600 


5 


0.285 


0.619 


8 


0.102 


L/L Edd 




H/3 


-0.036 


7 


0.939 


0.200 


9 


0.606 






Ha 


-0.300 


5 


0.624 


0.024 


8 


0.955 




FWHM 


H/3 


-0.286 


7 


0.534 


0.067 


9 


0.865 






Ha 


-0.200 


5 


0.747 


-0.071 


8 


0.867 


C iv Blueshift 


<ri 


H/3 


0.536 


7 


0.215 


0.033 


9 


0.932 






Ha 


0.100 


5 


0.873 


-0.214 


8 


0.610 




FWHM 


H/3 


0.679 


7 


0.094 


0.133 


9 


0.732 






Ha 


0.600 


5 


0.285 


0.143 


8 


0.736 


C IV assymetry 




H/3 


0.429 


7 


0.337 


0.117 


9 


0.765 






Ha 


0.300 


5 


0.624 


-0.048 


8 


0.911 




FWHM 


H/3 


-0.393 


7 


0.383 


-0.333 


9 


0.381 






Ha 


0.300 


5 


0.624 


-0.476 


8 


0.233 


AL A (1350A) / AL A (5100A) 


<ri 


H/3 


0.929 


7 


0.003 


0.883 


9 


0.002 






Ha 


1.000 


5 


0.000 


0.809 


8 


0.015 




FWHM 


H/3 


0.750 


7 


0.052 


0.767 


9 


0.016 






Ha 


0.700 


5 


0.188 


0.857 


8 


0.007 



Note. — The table shows the correlation strength of the BH mass residuals as a function of each 
different variable, quantified by the Spearman rank order coefficient, r s . Results are shown for group 
I and the combination of groups I and II measurements. In each case, N indicates the number of 
QSOs used to estimate the correlation strength and P ran indicates the probability of observing such a 
correlation by chance if the variables are uncorrelatcd. 
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Table 7. Linear Fits to Correlations of BH Mass Residuals with AGN color 



/~1 ft r 

Kj IV 


Balmer 


Group I Estimates 


Group I & II Estimates 


Width 


Line 


a 


b 




a 




b 


Prescription A 


















H/3 


0.64 ± 0.13 


-0.13 ± 0.06 


0.68 


± 


0.16 


-0.13 ± 0.08 




Ho 


0.58 ± 0.15 


-0.11 ± 0.06 


0.57 


± 


0.11 


-0.17 ± 0.06 


FWHM 


H/3 


0.89 ± 0.25 


-0.20 ± 0.12 


0.95 


± 


0.22 


-0.19 ± 0.11 




Ha 


0.75 ± 0.30 


-0.20 ± 0.13 


0.79 


± 


0.15 


-0.23 ± 0.08 


Prescription B 


















H/3 


0.60 ± 0.11 


-0.18 ± 0.05 


0.68 


± 


0.17 


-0.16 ± 0.08 




Ha 


0.51 ± 0.14 


-0.14 ± 0.06 


0.58 


± 


0.16 


-0.22 ± 0.08 


FWHM 


H/3 


0.86 ± 0.25 


-0.23 ± 0.12 


0.91 


± 


0.22 


-0.22 ± 0.10 




Ha 


0.72 ± 0.30 


-0.23 ± 0.13 


0.76 




0.16 


-0.27 ± 0.08 


Combined Sample 
















FWHM 


H/3 


0.82 ± 0.18 


-0.40 ± 0.07 


0.85± 0.18 


-0.40± 0.07 



Note. — The fits discussed in §4.21 correspond those performed using the prescription B 
C iv line-widths. Fits obtained using prescription A measurements are shown for complete- 
ness. The fits to the combined sample discussed in f|5]are also reported here. 
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Table 8. Correlations of Line- Width Estimates 



C IV Prescription 


Civ 


Balmer 


Group I Estimates 


Group I & II Estimates 




Width 


Line 


1 s 


N 


1 ran 


' s 


/V 


P 

1 ran 


Without color correction 


















Prescription A 


07 


H/5 


0.214 


7 


0.644 


0.300 


9 


0.433 






Ha 


0.700 


5 


0.188 


0.786 


8 


0.021 




FWHM 


H/5 


-0.143 


7 


0.760 


0.000 


9 


1.000 






Ha 


0.300 


5 


0.624 


0.500 


8 


0.207 


Prescription B 


07 


H/5 


0.679 


7 


0.094 


0.333 


9 


0.381 






Ha 


0.900 


5 


0.037 


0.524 


8 


0.183 




FWHM 


H/5 


-0.143 


7 


0.760 


0.000 


9 


1.000 






Ha 


0.300 


5 


0.624 


0.500 


8 


0.207 


Combined Sample 


FWHM 


H/5 


0.346 


52 


0.012 


0.342 


54 


0.011 


Combined Sample (S/N > 10) 


FWHM 


H/5 




31 


i 3 * in -3 


D 539 


33 


i 9 x m~ 3 


With color correction 


















Prescription A 


07 


H/5 


0.571 


7 


0.180 


0.550 


9 


0.125 






Ha 


0.900 


5 


0.037 


0.905 


8 


0.002 




FWHM 


H/5 


0.357 


7 


0.432 


0.383 


9 


0.308 






Ha 


0.400 


5 


0.505 


0.738 


8 


0.037 


Prescription B 


07 


H/5 


0.929 


7 


0.003 


0.467 


9 


0.205 






Ha 


0.900 


5 


0.037 


0.452 


8 


0.260 




FWHM 


H/5 


0.357 


7 


0.432 


0.383 


9 


0.308 






Ha 


0.400 


5 


0.505 


0.667 


8 


0.071 


Combined Sample 


FWHM 


H/5 


0.326 


52 


0.018 


0.323 


54 


0.017 


Combined Sample (S/N > 10) 


FWHM 


H/5 


0.602 


31 


3.4 x 10~ 4 


0.587 


33 


3.3 x 10- 4 



Note. — The table shows the correlation strength of the C iv and Balmer line- width estimates for both 
Civ line- width measurement prescription, quantified by the Spearman rank order coefficient, r s . Results 
are shown for group I and the combination of groups I and II measurements. In each case, indicates the 
number of lenses used to estimate the correlation strength and P TaD indicates the probability of observing 
such a correlation by chance if the variables are uncorrelated. 
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Fig. 1.— (Top) LUCIFER H- and i^-band spectra of SDSS1138+0314. The black solid line 
shows the spectrum obtained by performing the sky subtraction with the median combination 
of the sky frames while the gray line shows that obtained by using the modified version of 
the COSMOS software described in the text. (Bottom Left) Spectral region around H/3. 
Overlaid on top are the best fit continuum (black dotted line) and narrow (black short-dashed 
line) and broad line components (black long-dashed line), as well as their sum (black solid 
line) and the error spectrum (thin gray solid line). (Bottom Right) Same as bottom left but 
for Ha. 
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Fig. 2.— LUCIFER J-band spectrum of HS0810+2554 (gray solid line). Overlaid are the 
best fit continuum and Fell emission (black dotted line), narrow line emission (black short- 
dashed line), broad line component (black long-dashed line) and their sum (black solid line), 
as well as the error spectrum (thin gray solid line). 
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Fig. 3. — LIR I S nea r-IR spectra of images A and B of SBS0909+532 obtained by 



Mediavilla et al.l (120 lOh . The top panel shows the complete spectrum while the bottom 
four panels show the spectral regions around Ha and H/3 of each quasar image. Overlaid 
on top are the best fit continuum and narrow- and broad-line components, as well as their 
overall sum, using the same line-styles as in Figure [TJ 
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Fig. 4. — UV/optical spectra of all objects used in this study. For each object the panel 
shows the full spectrum, the source from which it was obtained and its redshift. 
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Fig. 4. — Continued 
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1500 1650 1500 1650 

Rest-frame A (A) 

Fig. 5. — The Figure shows the region of each UV/optical spectrum around Civ along with 
the best fit continuum (dotted line) around the line, the best fit line profile (long dashed) and 
the addition of both (short dashed). For each spectrum, the fits obtained using prescription 
A are shown in the left panel while those obtained following prescription B are shown in the 
right panel. We could not obtain good fits for SDSS1138+0314 and SBS0909+532 ClV lines 
(see Appendix lAl for details). 
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Fig. 5. — Continued 
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Fig. 5. — Continued 
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Fig. 6. — Comparison of C IV BH masses deri ved fro m the FWHM and ai velocity widths for 
all objects in our sample using the relations of lVPOa The left (right) panel compares the BH 
mass estimates based on the prescription A (B) line-width measurements of C IV. Masses are 
equal along the dashed line and the dotted line correspond to the best fit offset. The objects 
with arrows correspond to those for which we believe our C IV line- width measurements to 
be lower bounds. 
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Fig. 7. — Comparison between BH masses estimated from the prescription B o\ and FWHM 
of Civ and from the FWHM of Ha and H/3. For the estimates based on the line dispersion 
of Civ we have added the systematic offset of 0.24 dex described in §4.11 Solid symbols 
correspond to the objects with group I Ha or H/3 line-width estimates, while open sym- 
bols correspond to those wi th group II estimates. Six-pointed stars mark the objects not 
considered in the analysis of iGPLlOl . The dotted line shows where the BH masses are equal. 
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Fig. 8. — Ratio between the Civ and H/3/Ha mass estimates as a function of the corre- 
sponding hydrogen line mass estimate. Symbols and lines have the same definitions as in 
Figure 
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Fig. 9. — Same as Figure [5J but as a function of the UV continuum luminosity. 
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Fig. 10. — Same as Figure El but as a function of the 5100A continuum luminosity. 
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Fig. 11. — Same as Figure El but as a function of redshift. 
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Fig. 12. — Same as Figure El but as a function of the estimated Eddington ratio. We used a 
fa ctor of 11.91 to con vert between ALa(5100A) and LboI, as determined from the AGN SED 
of lAssef et all feoioh . 
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Fig. 13- 



Same as Figure El but as a function of the blueshift of the C IV line. 
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Fig. 14. — Same as FigureEl but as a function of the asymmetry of the C IV line, parametrized 
as the ratio of the widths red and blue of the line centroid. We do not show lower bounds 
on the C IV line- width as those object do not have a well defined blue side width due to the 
presence of absorption. 
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Fig. 15. — Same as Figure El but as a function of the ratio of the UV to optical continuum 
luminosities. The solid line shows the best fit linear relation to all objects with group I 
Balmer-line width estimates (solid symbols) and the dashed line shows the linear relation 
obtained when also including object with group II estimates (open symbols). Note that we 
do not include objects with lower bound C IV line widths on the fits. 
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Fig. 16. — Same as Figure [7J but after correcting the Civ BH masses for the dependence on 
the ratio of the UV to optical continuum luminosities observed in Figure [151 
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Fig. 17. — Comparison between the Civ and Balmer lines measured line-widths. Points and 
lines have the same meaning as in Figure [7J We do not show lower bounds for clarity. 
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Fig. 18. — Comparison between the ClV and Balmer lines measured line-widths after ap- 
plying correction from equation (jUJ). Points and lines have the same meaning as in Figure [7J 
We do not show lower bounds for clarity. 
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Fig. 19. — Panel a) sho ws the BH masses esti mated using eqns. (HI), (jSJ) and (Ej) for the 



and Dietrich et al. 



obj ects in the samples o f IVP06I (open squares), iNetzer et al.l ( 120071 ) (solid gray pentagons) 
J2009h ( s) for which this was possible (see $5]for details). 



Objects in our sample are shown by the solid and open six-pointed stars and circles, keeping 
the point style conventions used in previous plots. The dotted line shows where the masses 
are equal. Error-bars are not shown in order to make the plot more legible. Panel b) shows 
the results after applying the continuum slope correction from Table [7J 
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Fig. 20. — Residuals between BH masses estimated from the FWHM of the Civ and H/3 
broad emission lines as a function of t he log arithm of ratio of the continu u m lum inosities at 
1350A and 510 0A for the samples o f IVPOq {open squares), INetzer et al.l ( 20071 ) {solid gray 
pentagons) and iDietrich et al.l (120091 ) {solid black triangles) as well as our sample {solid and 
open six-pointed stars and circles). Point-styles have the same definitions as in Figure [T9l 
The solid line shows the best-fit linear relation for our data while the dashed line shows the 
best-fit to the combined sample. Error-bars are not shown in order to make the plot more 
legible. 
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Fig. 21. — Residuals between BH masses estimated from the FWHM of the Civ and B./3 
broad emission lines as a function of the logarithm of H/3-based BH mass (top left), redshift 
(top right), Eddington ratio (bott om left ) and continuu m luminosity at 5100 A( bottom right). 
Objects bel ong to the samples o f IVP06I (open squares), iNetzer et al. I 653) ( solid gray pen- 
tagons) and lDietrich et al.l ( 120091 ) (solid black triangles) as well as our sample (solid and open 
six-pointed stars and circles). Point-styles have the same definitions and in Figure [T9"l The 
dotted line shows where the masses are equal. Error-bars are not shown in order to make 
the plot more legible. 
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Fig. 22. — Comparison of the Civ and H/3 line-widths for the sample we have compiled 
from the literature. Point types and line styles are the same as for Fig. [191 The large 
gray hexagon shows SDSS1151+0340. Top panels show the complete literature sample while 
bottom panels only show objects with spectra that have continuum S/N > 10 in the vicinity 
of Civ. 
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Fig. 23.— Spectra of the QSO SDSS1 151+0340 obtained by SDSS (top), with MDM/OSMOS 
(middle) and with Palomar /Double Spectrograph^ bottom). The arrows mark the probable 
absorption troughs near the Civ line. The spectra have been resampled to a common 
resolution, and the continuum S/N shown in the upper left corner of each panel has been 
calculated in the same way as for all other objects in our sample. 



